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POLARIZATION  MULTIPLEXING  SYSTEMS 


Thomas  D.  Moran  and  Bernard  J.  Lamberty 
Boeing  Electronics  Company,  High  Technology  Center 


Abstract 

In  a  polarization  diversity  system  it  is  conceptually  feasible  to  multiplex  two 
orthogonally  polarized  signals  onto  a  single  channel.  The  savings  in  system 
hardware  could  be  significant,  particularly  if  the  antenna  is  a  phased  array  and 
multiplexing  can  be  achieved  near  the  front  end  of  each  element.  Multiplexing 
however,  increases  T/R  module  complexity  and  could  degrade  polarization 
orthogonality  below  acceptable  levels  if  isolation  between  channels  is  not 
preserved.  A  study  was  conducted  to  evaluate  advantages  and  disadvantages  of 
three  different  multiplexing  techniques;  frequency,  quadrature,  and 
pseudorandom  code.  Increase  in  bandwidth,  phase  sensitivity,  and  complexity 
relative  to  a  two  channel  system  are  compared. 

1.0  INTRODUCTION 

A  polarization  diversity  system  receives  and  processes  two  orthogonally 
polarized  signals  simultaneously.  This  normally  requires  two  independent  sets  of 
hardware  throughout  the  system  up  to  the  location  where  the  two  signals  are 
processed  or  combined.  Depending  on  the  specific  application,  a  minimum  level 
of  isolation  along  with  phase  and  amplitude  balance  must  be  preserved  between 
the  two  signals  to  maintain  polarization  orthogonality.  If  these  two  signals  could 
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be  multiplexed  onto  a  single  channel  while  maintaining  adequate  channel  isolation 
and  balance,  a  significant  system  simplification  and  therefore  cost  savings  is 
potentially  available. 

Radar,  communications,  and  ECM  are  examples  of  systems  which  often  use 
polarization  diver sity  and  so  could  benefit  from  polarization  multiplexing.  An  L- 
band  radar  application  was  used  to  evaluate  this  concept.  To  provide  parameter 
values  for  comparison  between  multiplexing  candidates,  the  following 
specifications  were  selected  for  this  study:  channel  isolation  >30  dB,  frequency 
band  center  =  1  GHz,  bandwidths  =  10  MHz  and  100  MHz.  Extrapolations  to 
other  applications  can  be  made  based  on  results  of  this  study. 

Figure  1  shows  schematic  diagrams  of  typical  T/R  modules  used  to  evaluate 
the  polarization  multiplexer  approaches.  Figure  la  shows  a  module  without 
polarization  multiplexing;  Figure  1b  shows  the  same  module  with  multiplexing 
incorporated.  Note  that  transmission  is  on  a  single  polarization  whereas  the 
return  signal  is  received  simultaneously  on  orthogonally  polarized  channels. 
Such  a  system  might  be  used  to  suppress  interfering  signals  via  polarization 
nulling  or  to  characterize  target  scattering  as  a  function  of  polarization.  The  dual 
channel  system  requires  at  least  two  combiner  networks  compared  to  one  for  a 
multiplexed  system.  The  potential  to  eliminate  one  array  combiner  and  its  phase 
matched  cables  seemed  to  outweigh  the  complexity  added  by  a 
multiplexer/demultiplexer  system  -  especially  if  the  multiplexer  could  be 
implemented  with  a  microwave  monolithic  integrated  circuit  (MMIC). 
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2.0  CANDIDATES 


The  three  multiplexing  candidates  are  frequency,  quadrature  and 
pseudorandom  code.  These  candidates,  shown  in  Figure  2  along  with  their 
respective  frequency  spectra  in  Figure  3,  are  described  below. 

In  frequency  division  multiplexing  a  common  local  oscillator  is  mixed  with 
each  polarization  channel.1  The  upper  sideband  (fRF  +  fLO)  is  used  by  ,he 
vertically  polarized  channel  and  the  lower  sideband  (fRF  -  fLO)  by  tbe 

horizontally  polarized  channel.  These  two  frequency  bands  are  summed  and  sent 
to  the  N-way  combiner  along  a  common  transmission  line.  Demultiplexing  is 

achieved  at  the  array  output  by  power  splitting  and  filtering  to  separate  the 
upper  and  lower  sidebands.  These  sidebands  are  then  mixed  with  fLO  and  filtered 

to  recover  the  orthogonally  polarized  signals. 

In  quadrature  division  multiplexing  the  horizontal  channel  is  mixed  with  the 
local  oscillator  and  the  vertical  channel  with  the  local  oscillator  phase  shifted  90 
degrees.2  This  places  the  mixer  products  (both  upper  and  lower  sidebands)  of 
channel  one  and  two  in  phase  quadrature  allowing  both  channels  to  be  combined 
onto  one  line  while  preserving  isolation.  Demultiplexing  is  accomplished  by 
power  splitting  the  common  line  and  re-mixing  with  the  in-phase  and  quadrature 
local  oscillator  to  preserve  the  time  alignment  of  each  channel  at  the 
multiplexer.  The  mixer  products  are  such  that  cancellation  occurs  eliminating 
the  undesired  channel  leaving  the  desired  signal  at  the  mixer  output.  Bandpass 
filtering  is  used  to  remove  extraneous  products  incurred  through  the 
demultiplexing  process. 
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In  code  division  multiplexing  each  channel  is  multiplied  by  a  unique 
pseudorandom  code  before  the  channels  are  summed  and  sent  along  a  single 
transmission  line.^-4  Typically  maximum  linear  code  sequences  are  chosen  to 
encode  the  channels.  The  same  code  sequence  can  be  applied  to  each  mixer  within 
the  multiplexer  provided  that  one  channel's  code  is  delayed  by  at  least  one  code 
bit.  Since  code  sequences  of  this  type  have  a  DC  component,  isolation  is  limited 
by  code  length.  This  DC  component  is  inversely  proportional  to  code  length  so  if 
sufficiently  long  codes  are  chosen  a  low  DC  component  results.  For  example,  a 
31 -bit  code  has  a  sufficiently  small  DC  component  to  allow  nearly  30  dB  of 
isolation  between  channels.  Since  the  code  is  pseudorandom  the  spectral 
distribution  is  sin(x)/x  with  main  lobe  nulls  at  plus  and  minus  the  code  rate 
about  the  center  frequency.  For  a  10  MHz  information  bandwidth  a  310  MHz  code 
rate  is  needed  resulting  in  a  620  MHz  multiplexed  channel  bandwidth. 
Demultiplexing  is  accomplished  by  mixing  the  codes  with  the  multiplexed  channel 
in  such  a  way  as  to  maintain  the  same  time  relationship  that  occurred  at  the 
multiplexer.  A  bandpass  filter  is  then  used  to  separate  the  desired  channel  from 
the  sin(x)/x  spectrum  of  the  undesired  channel. 


2.1  IMPLEMENTATION  DETAILS 

Single  sideband  frequency  division  multiplexing  can  be  implemented  within 
each  T/R  module  in  two  basic  ways,  through  use  of  either  single  sideband 
modulators  or  single  sideband  filters.  The  single  sideband  modulator  approach 
requires  four  mixers  and  two  0-180  degree  hybrids  per  module.  The  single 
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sideband  filter  approach  is  more  direct  but  requires  phase  matched,  flat  group 
delay  filters  within  each  T/R  module. 

Quadrature  division  multiplexing  is  simpler  from  a  hardware  standpoint, 
although  to  achieve  the  desired  channel  isolation,  phase  match  between 
multiplexer  and  demultiplexer  local  oscillators  must  be  tightly  controlled  (±1 .8 
degrees  for  30  dB  channel  isolation).  Additionally,  this  local  oscillator  phase 
sensitivity  can  have  an  effect  on  selection  of  a  T/R  module  configuration.  The 
common  transmit/receive  phase  shifter  shown  in  Figure  1b  alters  this  phase 
match  by  delaying  the  multiplexed  channel  with  respect  to  the  local  oscillator 
thus  degrading  isolation.  For  L-band  systems  of  narrow  bandwidth  (=1%)  this 
effect  can  be  tolerated  because  the  local  oscillator  frequency  can  be  very  low.  For 
wideband  systems  (=10%)  however,  isolation  degradation  is  sufficiently  severe 
that  if  a  phase  shifter  is  to  be  included  within  the  multiplexed  channel  (as  shown 
in  Figure  1b),  phase  compensation  in  the  demultiplexed  local  oscillator  is 
required.  Alternately,  two  receiver  phase  shifters  could  be  placed  ahead  of  the 
multiplexer  such  that  the  multiplexed  channel  and  local  oscillator  remain  in  time 
synchronization  regardless  of  phase  shifter  position.  This  would  require  an 
additional  phase  shifter  and  thus  would  increase  array  complexity. 

Code  division  multiplexing  is  similar  to  the  quadrature  technique  in  terms  of 
T/R  module  complexity.  While  this  technique  requires  a  phase  matched  local 
oscillator  like  the  quadrature  case,  the  isolation  is  much  more  tolerant  of  phase 
mismatch.  However  a  much  higher  frequency  "local"  oscillator  must  be  used 
which  offsets  this.  A  potential  advantage  of  this  method  is  that  more  than  two 
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channels  may  be  multiplexed.  This  could  be  used  to  advantage  in  a  dual 
polarization  system  using  sum  and  difference  channels  for  angle-of-arrival  data. 

2.2  MEASUREMENTS 

Measurements  were  made  of  two  of  the  multiplexed  approaches,  quadrature  and 
code,  to  verify  that  the  specified  level  of  isolation  was  achievable.  No 
measurements  were  made  of  the  frequency  technique  because  its  performance  is 
easily  predicted. 

A  model  of  a  quadrature  multiplexed  system  was  measured  with  the  test  setup 
shown  in  Figure  4a  and  pictured  in  Figure  5.  Slightly  different  frequencies  were 
used  fcr  each  channel  so  that  relative  channel  isolation  could  be  displayed  on  a 
spectrum  analyzer.  Measurements  were  made  with  a  51  MHz  local  oscillator  and 
1000  MHz  and  1003  MHz  carriers  for  each  of  the  two  polarization  channels.  The 
51  MHz  local  oscillator  was  used  to  accommodate  the  available  5%  bandwidth  1 
GHz  tunable  bandpass  filters.  (Note:  isolators  in  the  demultiplexer  are 
necessary  to  prevent  high  order  products  reflected  form  the  bandpass  filter  from 
reducing  channel  isolation  via  recombination  within  the  mixer.)  Measurements 
indicated  that  30  dB  of  isolation  could  easily  be  obtained  with  the  quadrature 
multiplexer  method. 

The  test  setup  for  the  code  division  multiplexer  model  is  similar  to  that  of  the 
quadrature  setup  except  the  0-90  degree  hybrids  are  not  required  and  no  output 
filtering  was  attempted  due  to  the  low  code  rate  necessitated  by  available 
equipment.  A  [5,2]  31 -bit  maximum  linear  code  with  a  10  MHz  chip  (clock) 
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rate  was  used  to  drive  the  mixers  as  shown  in  Figure  4b.  A  digital  delay  circuit 
was  used  to  generate  a  one-bit  delay  for  channel  2.  Isolation  measured  =30  dB  as 
predicted  for  the  31 -bit  code. 

2.3  COMPARISON  OF  ADVANTAGES  AND  DISADVANTAGES 

The  key  advantage  of  a  frequency  division  multiplexed  system  is  complete 
insensitivity  to  multiplexer/demultiplexer  local  oscillator  phase.  Its  most 
important  drawback  is  that  accurate  phase  matched  filters  both  at  the 
multiplexer  and  demultiplexer  may  prove  difficult  to  implement. 

The  key  advantages  of  the  quadrature  division  multiplexer  approach  are  that 
the  multiplexer  within  the  T/R  module  is  simple  and  easy  to  phase  match.  The 
low  frequency  5  MHz  local  oscillator  and  its  quadrature  can  be  generated  digitally 
with  a  clock.  Phase  matching,  while  critical  is  simplified  due  to  the  low 
frequencies  involved.  This  technique  tends  to  work  best  for  narrow  bandwidth 
systems  although  it  can  be  adapted  for  wide  bandwidth  operation  by  adding 
additional  hardware. 

Code  division  multiplexing  is  much  less  phase  sensitive  than  quadrature 
multiplexing.  However,  because  the  code  rate  for  a  typical  system  with  1% 
bandwidth  is  more  than  50  times  that  of  the  quadratuie  system,  the  low  phase 
sensitivity  of  this  approach  is  overshadowed  by  the  other  technique's  lower  local 
oscillator  frequency.  As  a  result  of  the  high  local  oscillator  frequency,  the 
multiplexed  channel  bandwidth  is  quite  high  -  nearly  1  GHz  for  a  1%  bandwidth 
system  and  much  higher  for  a  10%  bandwidth  system.  Figure  6  compares 


isolation  vs.  local  oscillator  line  length  variation  for  the  quadrature  and  code 
division  techniques. 

Figure  7  summarizes  advantages  and  disadvantages  of  the  three  multiplexing 
alternatives. 

3.0  CONCLUSIONS 

In  a  two  channel  radar  a  multiplexed  system  has  the  potential  for  reducing  the 
number  of  high  frequency  power  combiners  and  associated  phase  matched  cables 
by  a  factor  fo  two.  Thi*  is  achieved  at  the  expense  of  a  more  complex  T/R  module 
and  the  distribution  of  phase  matched  local  oscillator  lines. 

For  the  L-band  radar  example  with  1%  bandwidth  the  quadrature 
multiplexer  approach  appeared  optimal.  Its  low  frequency  local  oscillator 
requirement,  low  channel  bandwidth,  and  simple  module  multiplexer  held  the 
greatest  promise  despite  its  phase  sensitivity.  Monolithic  active  mixers  have  the 
potential  for  yielding  small,  inexpensive,  phase  matched  multiplexers,  which 
because  of  their  MMIC  implementation  could  improve  the  cost/pertormance  ratio 
of  a  phased  array  system. 
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Figure  3.  Spectral  Distribution  for  Dual  Polarization 
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Quadrature  multiplexer 
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Figure  4.  Measurement  Test  Setup 
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Figure  5.  Quadrature  Division  Multiplexer/Demultiplexer  Test  Setup 
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INTERFERENCE  SUPPRESSION  USING  AN  ADAPTIVE  POLARIZATION  COMBINER 


Bernard  J.  Lamberty  and  Robert  P.  Friedman 
Boeing  Electronics  Company,  High  Technology  Center 


ABSTRACT 

Polarization  diversity  systems  receive  and  process  dual,  orthogonally 
polarized  signals.  A  circuit  consisting  of  hybrids  and  phase  shifters  can  be  used 
for  combining  these  signals  which  may  be  of  arbitrary  magnitude  and  phase. 
With  proper  adjustment  of  the  phase  shifters,  either  open-loop  or  adaptively, 
total  input  power  is  sent  to  one  output  port  while  a  signal  null  occurs  at  the 
other.  Thus  the  circuit  is  suited  for  either  optimizing  received  power  from  a 
desired  signal  or  for  suppressing  undesirable  in-band  interfering  signals.  A 
study  was  conducted  to  evaluate  this  polarization  combining  network  including 
determination  of  effects  of  component  errors  on  null  depth.  Of  particular 
interest  were  quantization  effects  of  using  digital  phase  shifters  and  amplitude 
detectors  in  the  circuit.  Null  depth  as  a  function  of  amplitude  and  phase  of  input 
signals  was  calculated  and  compared  with  measurements  for  combiners  using 
both  4-  and  5-bit  phase  shifters  and  for  amplitude  detection  quantization  levels 
of  up  to  3  dB.  Results  are  used  to  predict  effectiveness  of  interfering  signal 
suppression  as  a  function  of  incident  signal  polarization. 
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1.0  INTRODUCTION 


Polarization  diversity  can  be  used  for  signal  enhancement  or  to  suppress  an 
interfering  signal  whose  polarization  is  different  from  the  desired  signal.  Figure 
1  is  a  schematic  of  a  simple  network  that  can  be  used  to  do  this.  V1  and  V2  are 

the  incident  signal  voltages  received  by  the  orthogonally  polarized  antennas  of  the 
diversity  system.  These  voltages  are  of  arbitrary  relative  magnitude  and  phase. 
The  first  phase  shifter,  PH1 1 ,  is  used  to  equalize  their  phases.  When  these  equal 
phase  voltages  are  passed  through  the  first  hybrid,  HI ,  its  output  voltages  are  of 
equal  magnitude.  The  second  phase  shifter,  PHI2,  is  used  to  adjust  these  equal 
magnitude  voltages  to  be  orthogonal  in  phase.  Finally,  when  the  equal  magnitude, 
orthogonal  phase  voltages  pass  through  the  second  hybrid,  H2,  they  add  in  phase 
at  one  output  and  cancel  at  the  other.  Thus,  the  combiner  can  be  used  to  either 
extract  all  the  power  from  an  arbitrarily  polarized  signal  at  the  sum  output 
port,  or  to  cancel  an  in-band  interfering  signal  at  the  null  output  port.*  This 
paper  concentrates  on  the  level  of  cancellation  available  from  practical 
polarization  combiner  circuits. 

2.0  DETAILS  OF  OPERATION 

Voltage  magnitude  and  phase  relationships  at  successive  junctions  of  the 
polarization  combiner  circuit  shown  in  Figure  1  are: 

V1  =V1eie°,  V^-V^00,  V2  =  V2ei0°  (1) 

For  processing,  the  interfering  signal  must  have  some  parameter  to  distinguish 
it  from  the  desired  signal  (e.g.  bandwidth) 
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(2) 


v1  =  1  (Viei0°4.  Vce'iTC/2^ 

V2" 

(I 

V2  =  _J_  (V1e')7t/2+  V2ei°°)  (3) 

< 2 

ill 

V1  =  _J_  (V1  e')“+  V2e'i(a  +  7t/20  (4) 

VpQ  =  _1_  0^6*1“+  V2e'j<a  +  7t/2)  +  V1e'j;c+  V2e‘iIt/2h  (5) 

VgQ  =  __i_  (V1e'j(a+It/2)+  V2e'Ha  +  7t)  +  v1e'H7t/2)+  V2e'i°°)  (6) 

4 2 


To  find  the  values  of  alpha  that  sends  all  the  input  power  to  Port  A  and  zero  power 
to  Port  B,  the  real  and  imaginary  parts  of  VBq  are  each  set  equal  to  zero.  This 


(7) 


(8) 
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Alpha,  plotted  as  a  function  of  lv1/V2l  in  Figure  2,  varies  from  0  degrees  for 
V i  =0  to  90  degrees  for  V.|=V2  to  180  degrees  for  V2=0.  If  the  values  of  alpha 

shown  in  Figure  2  are  increased  by  180  degrees,  maximum  power  appears  at 
VgQ  and  the  null  appears  at  V^q-  Relative  phase  between  input  signals  is 

arbitrary  so  PHI1  must  provide  360  degrees  of  phase  shift. 


3.0  ADAPTIVE  COMBINER  EVALUATION 

If  the  polarization  combiner  is  to  be  adjusted  or  controlled  remotely,  digital 
phase  shifters  may  be  more  convenient  for  PHI1  and  PHI2  than  continuously 
variable  designs.  Then  since  only  a  limited  number  of  discrete  values  of  phase 
shift  are  available,  total  signal  cancellation  can  only  be  achieved  for  certain 
specific  input  signal  amplitude  ratios  and  relative  phases.  (See  Figure  2).  The 
null  depth  for  other  input  signal  combinations  is  a  function  of  number  of  bits  in 
the  phase  shifters. 

In  addition,  if  the  combiner  is  adaptive,  the  phase  shifters  must  be  driven  by 

closed  loop  circuits  as  shown  conceptually  in  Figure  3.  Here  the  phase  detector 
circuit  drives  the  first  phase  shifter,  PHI1,  until  V1  and  V2  are  in-phase  or  as 

close  to  in-phase  as  that  digital  phase  shifter  will  allow.  Then  the  amplitude 

detection  circuit  drives  the  second  phase  shifter,  PHI2,  until  a  minimum  signal 
is  present  at  Vg.  In  another  concept,  the  amplitude  detection  circuit  is  located 

between  PHI1  and  HI.  PHI2  is  then  driven  to  its  appropriate  phase  setting  based 
on  stored  lookup  tables  taken  from  Figure  2.  Because  these  or  other 
implementations  could  involve  measurement  of  lv1/V2l,  effects  of  amplitude 

measurement  resolution  on  null  depth  were  examined  for  quantization  levels  of  1 
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dB,  2  dB  and  3  dB.  Hybrid  and  phase  shifter  error  effects  can  also  affect  null 
dep  h  and  so  were  included  in  the  study. 

3.1  PASSIVE  COMPONENTS  ERROR  EFFECTS  ANALYSIS 

The  sources  of  error  in  the  adaptive  polarization  combiner  circuit  are;  (1) 
phase  shifter  errors  (quantization,  phase  detection,  and  other  deviations  from 
ideal  phase  settings),  (2)  input  voltage  ratio  amplitude  detection  errors,  (3)  90 
degree  hybrid  amplitude  and  phase  imbalances,  (4)  phase  shifter  insertion  loss 
variations,  (5)  phase  shifter  and  hybrid  VSWR  and  (6)  hybrid  isolation.  These 
error  terms  are  listed  in  decreasing  order  of  importance  under  the  assumption 
that  the  combiner  uses  typical  hybrids  and  digital  phase  shifters  with  a  least 
significant  bit  of  a  few  degrees.  Equations  1-8,  expanded  to  include  the  first  4  of 
these  error  sources,  were  used  to  predict  null  depths  for  input  voltage  ratios  of  0 
to  60  dB.  Effects  of  the  most  important  sources  were  calculated  first  with  less 
important  sources  added  in  successive  iterations. 

Figure  4  shows  null  depth  (dB  below  toial  input  power)  as  a  function  of  input 
voltage  ratio  if  V1  and  V2  are  in  phase.  Note  that  PHI2  is  the  only  source  of 

quantization  phase  error  for  this  case  since  PHI1  phase  =  0.  There  are  four 
shallow  nulls  of  approximately  20  dB  each  when  a  4-bit  phase  shifter  is  used 
(22.5  degree  phase  increment)  and  eight  shallow  nulls  of  approximately  26  dB 
each  when  a  5-bit  phase  shifter  is  used  (11.25  degree  phase  increments).  The 
deep  nulls  between  each  of  these  correspond  to  input  voltage  magnitudes  where 
the  required  phase  for  optimum  nulling  is  exactly  equal  to  one  of  the  phase 
shifter  bit  settings.  Shallowest  nulls  occur  when  the  required  phase  for  optimum 
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nulling  is  half  a  bit  from  the  two  nearest  available  phase  shifter  bit  settings. 
Shallow  nulls  are  lower  by  6  dB  for  every  additional  bit  in  the  phase  shifter. 
Identically  shaped  curves  exist  for  input  voiiage  magnitude  ratios  of  0  to  -60  dB. 
(See  Figure  2). 

Figure  5  is  the  same  as  the  4-bit  case  in  Figure  4  except  that  the  input 
voltage  phase  differences  range  from  0  to  11.25  degrees.  PH1 1 ,  used  to  co-phase 
these  voltages,  now  can  introduce  an  additional  quantization  error  which  is 
greatest  at  input  voltage  phase  differences  halfway  between  phase  shifter  bits 
(11.25  degrees).  Here  the  first  deep  null  degrades  to  20  dB  and  the  first  shallow 
null  to  17  dB.  Degradations  of  both  the  shallow  and  deep  nulls  are  less  severe  as 
input  voltage  ratios  increase.  However,  all  nulls  are  less  than  30  dB.  Similar 
curves  were  computed  for  the  5-bit  case.  Nulls  were  deeper  by  about  6  dB. 
Responses  are  periodic  with  input  signal  phase  difference.  The  periodicity  is 
related  to  quantization  bit  increments,  that  is,  the  same  response  occurs  for 
input  voltage  phase  differences  of  Nx(bit  increment)  with  N  an  integer. 

Figure  6  is  similar  to  the  case  shown  in  Figure  4  except  that  amplitude  ratio 

measurement  quantization  errors  are  included.  Here  the  decision  to  switch  phase 
shifter  settings  is  made  in  1  dB,  2  dB  or  3  dB  increments  of  I V  ^ /V2 1  -  The 

resulting  null  depth  is  compared  to  0  dB  amplitude  resolution.  Using  the  2  dB 
increment  as  an  example,  if  IV-j/N^I  is  between  4  dB  and  7  dB,  PHI2  is  set  as  if 

that  ratio  were  6  dB.  This  degrades  null  depth  for  I /V2 1  between  5  and  5.4 
dB,  since  the  phase  shifter  would  have  a  different  setting  in  that  range  for  zero 
amplitude  resolution  error.  For  the  4-bit  phase  shifter  case,  this  range  of 
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amplitude  resolution  errors  affects  the  shallow  nulls  only  which  are  degraded  by 
0.9  dB,  2.2  dB  and  3.7  dB  respectively. 

Similar  calculations  were  made  for  5-bit  phase  shifters  in  the  combiner. 
Shallow  nulls  degrade  by  about  3  dB  for  both  the  1  dB  and  2  dB  amplitude 
resolution  cases  and  by  over  6  dB  for  the  3  dB  resolution  case.  Also,  the  3  dB 
resolution  case  results  in  the  loss  of  two  deep  nulls  completely.  Further  null 
degradation  occurs  if  the  input  voltage  phases  are  separated  by  half  a  bit.  Then 
the  worse  case  null  is  about  16  dB  for  the  4-bit  case  and  about  19  dB  for  the  5- 
bit  case. 

Figure  7  shows  the  added  error  effects  of  hybrid  amplitude  imbalances  and 
phase  shifter  loss  variations.  Measured  values  of  individual  components  were 
used  for  these  predictions.  As  shown  in  the  figure,  results  depend  on  which  input 

voltage  is  larger  since  hybrid  imbalance  is  constant.  If  its  insertion  loss  is 
greater  on  the  V1  side,  the  effect  will  be  different  for  V^N/2  than  for  V^V 2  and 

vice  versa. 

3.2  MEASUREMENT  OF  ERROR  EFFECTS  IN  PASSIVE  COMPONENTS 

The  error  analyses  were  verified  by  measurements  using  the  polarization 
combiner  model  and  test  setup  shown  in  Figures  8  and  9.  Tests  were  conducted 
over  the  950  MHz  to  1150  MHz  band.  A  single  RF  source  feeds  a  3  dB  power 
divider  which  in  turn  feeds  the  sum  and  null  lines  in  the  test  setup.  In  the  sum 
line  there  is  a  discrete  level  variable  programmable  attenuator  with  attenuation 
bits  of  1,  2,  4,  10,  20  and  40  dB.  a  variable  analog  phase  shifter  is  used  lo  set 
input  signal  phase  differences.  Pads  are  used  between  components  to  minimize 
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VSWR  interactions.  The  model  uses  four  digital  phase  shifters,  but  only  the  two 
in  the  sum  line  are  used  to  adjust  phase.  The  two  difference  line  phase  shifters 
are  always  set  to  0  degrees  and  are  only  used  to  balance  insertion  loss  and  phase 
of  the  operated  phase  shifters. 

Output  null  depths  were  measured  from  950  MHz  to  1150  MHz  on  the  HP 
8573  network  analyzer  with  a  marker  on  1000  MHz  to  record  data  for  single 
frequency  analysis  comparison  verification. 

Figure  10  shows  a  comparison  of  measured  null  depths  versus  those  computed 
from  the  single  frequency  analyses.  This  plot  is  for  an  input  signal  phase 
difference  of  80  degrees  and  input  magnitude  ratios  of  0  to  60  dB.  All  measured 
component  error  terms  are  included  in  the  analysis.  The  agreement  is  excellent 
between  measured  and  computed  values.  Similar  comparisons  for  input  phase 
difference  of  0  and  40  degrees  also  showed  excellent  agreement. 

Figure  11  shows  a  comparison  of  measured  and  calculated  data  over  a  20% 
frequency  band  for  an  input  voltage  amplitude  ratio  of  10  dB  and  phase  difference 
of  40  degrees.  The  calculated  data  predicts  the  deepest  null  tu  occur  at  a 
frequency  greater  than  1000  MHz  which  is  in  agreement  with  the  measured  data. 

4.0  CONCLUSIONS 

A  simple  polarization  combiner  can  be  used  in  a  polarization  diversity  system 
to  suppress  interfering  signals  where  the  level  of  suppression  depends  on  the 
quality  of  the  components. 

Analysis  can  predict  combiner  performance  based  on  realistic  specifications 
and/or  measurements  of  its  components. 
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Digital  phase  shifters  and  amplitude  detection  circuits  may  be  required  to 
make  the  combiner  adaptive.  Over  17  dB  of  interference  suppression  can  be 
attained  using  4-bit  phase  shifters  and  over  23  dB  using  5-bit  phase  shifters. 

The  two  largest  contributors  to  degraded  nulls  are  the  quantization  errors  of 
these  digital  phase  shifters  and  amplitude  ratio  detectors. 
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Figure  1.  Polarization  Combiner  Network  and  Vector  Diagram 


Figure  2.  Phase  Shift  Required  for  Unequal  Signal  Combining 


Null  depth  -dB 


Figure  3.  Adaptive  Polarization  Combiner 
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Figure  4.  Polarization  Combiner  Predicted  Null  Depths;  4-Bit  and  5-Bit 
Phase  Shifter;  Input  Voltages  in  Phase 
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Figure  5.  Polarization  Combiner  Predicted  Null  Depth ;  4-Bit  Phase 
Shifters;  Input  Voltage  Phases  Variable 


|Vi/V2|-dB 

Figure  6.  Predicted  Null  Depth;  4-Bit  Phase  Shifters;  Input  Voltages  in 
Phase;  Amplitude  Detection  Levels  Variable 
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jure  10.  Comparison  Between  Calculated  and  Measured  Null  Depths; 
Combiner  With  4-Bit  Phase  Shifters;  80° Input  Voltage  Phase 
Difference;  1  GHz 


Figure  11.  Measured  and  Calculated  Null  Depths.  [V1/V2]  =  10  dB, 
at  -h  40° Phase 
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ABSTRACT 


PERFORMANCE  CHARACTERISTICS  OF  NOTCH  ARRAY 
ELEMENTS  OVER  A  6/1  FREQUENCY  BAND 

George  J.  Monser 

This  paper  presents  pattern  and  gain  characteristics  for  a 
notch  array  element  over  a  6  to  1  frequency  band. 

It  is  shown  that  the  lower  frequency  limit  occurs  indepen¬ 
dently  of  the  notch  design,  and  that  the  upper  frequency  limit 
depends  upon  the  element  design. 

Refinements  to  the  design  are  discussed,  leading  to  a  design 
operable  from  2.5  to  19.0  GHz. 


264 


PERFORMANCE  CHARACTERISTICS  OF  NOTCH  ARRAY 
ELEMENTS  OVER  A  6/1  FREQUENCY  BAND 


George  J.  Monser 

1.0  INTRODUCTION 

The  dual-polarized  array  consisting  of  printed-notches  in 
the  walls  of  the  horn  elements  was  reported  in  the  literature^ 
where  good  performance  characteristics  were  displayed  over  the 
continuous  frequency  band  from  4  to  18  GHz.  Since  the  printed 
notches  appeared  to  have  an  even  greater  bandwidth,  their  band¬ 
width  potential  was  checked. 

This  paper  shows  that  the  useful  bandwidth  for  this  type 
of  printed-notch  array  is  greater  than  6/1  (i.e.,  3  to  19  GHz); 
with  further  design,  performance  from  2.5  to  19  GHz  may  be 
realized . 

2.0  THE  TEST  MODEL 

Figure  1  shows  a  photograph  of  the  dual-polarized  array  with 
the  notches  built  into  the  walls  of  the  horn  elements.  Each 
consists  of  four  notches  fed  by  a  corporate  power-divider  struc¬ 
ture  as  shown  in  Figure  2.  The  transformers  within  the  dividers 


1.  Monser,  G.J.  (1984)  Considerations  for  Extending  the  Band- 
widths  of  Arrays  Beyond  Two  Octaves.  Paper  presented  at  the 
Antenna  Application  Symposium. 


were  designed  to  work  best  over  the  4  to  18  GHz  range  so  that 
part  of  the  element  rolloff  below  4  GHz  is  due  to  added  divider 
losses  which  is  further  aggravated  by  increased  coupling  between 
the  notch,  sub-array  elements. 

3.0  TEST  RESULTS 

Two  regions  of  bandwidth  extension  are  reported;  18  to 
19  GHz  and  3  to  4  GHz.  These  extensions  are  then  merged  with  the 
existing  band. 

3.1  ABOVE  18  GHz 

Figures  3  through  5  show  patterns  and  swept  gain  on  antenna 
boresight . 

The  array-plane  patterns  (smooth  and  broad)  (Figure  3)  shows 
that  adequate  coverage  is  provided  for  use  as  an  array  element. 

In  the  elevation  plane  the  patterns  (Figure  4)  are  similarly 
acceptable . 

Figure  5  shows  that  the  gain  is  a  smooth  function  with  no 
drop-outs  in  gain. 

3.2  BELOW  4  GHz 

Figures  6  through  9  show  patterns  and  swept  gain  on  antenna 
boresight . 

The  array-plane  patterns  (Figure  6)  show  that  adequate  broad 
coverage  is  provided  for  use  as  an  array  element. 


266 


Figure  7  shows  that  the  elevation  coverage  is  acceptable. 

Figure  8  shows  that  the  gain  is  a  smooth  function  with 
no  drop-outs  in  gain. 

Active  VSWR  for  this  element  is  shown  in  Figure  9  showing 
that  reasonably  good  transmission  efficiency  results  even  with 
the  electrically  close  spacing  between  array  elements. 

3 . 3  FULL  BAND 

Figure  10  shows  antenna  element  gain  over  the  full  band 
from  3  to  19  GHz. 

Figure  11  shows  the  elevation  coverage  ( orthogonal-to- 
array)  from  3  to  19  GHz. 

4.0  NOTCH  DESIGN 

The  basic  building  block  in  the  assembly  shown  in  Figure  2 
is  the  printed  circuit  notch.  Each  notch  is  excited  by  a  strip¬ 
line,  button-hook,  capacitively  coupled  line  (Figure  12).  The 
feed-line  cross-over  position  (distance  to  the  notch)  as  well 
as  the  extension  after  cross-over  and  geometry  (disk)  were 
empirically  determined  for  best  active  VSWR  corresponding  to 
a  scan  angle  of  _+  30°.  Both  active  and  passive  VSWR’s  were 
measured  over  the  4  to  18  GHz  band. 

The  slot-to-notch  transition  was  then  developed  from  a 
step  to  a  flare-step  configuration  (Figure  13)  which  showed 
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Improved  passive  match  without  affecting  the  active  match.  A 
total  of  four  notches,  spaced  0.42  inch  between  centers  and 
fed  by  a  four-way,  equal-split  divider  comprise  each  array 
element.  The  equal  power  division  was  selected  to  yield  the 
highest  gain  and  two  section  transformers  were  designed  and 
incorporated  in  the  corporate  feed  to  assure  a  good  match  at 
the  input. 

5.0  DISCUSSION  OF  RESULTS 

The  upper  frequency  limit  of  19  GHz  is  primarily  due  to 
pattern  degradation  although  the  transformers  in  the  corporate 
feed  are  beyond  their  design  band.  The  lower  frequency  limit  of 
3  GHz  is  primarily  due  to  poor  active  match  and  low  antenna  gain. 
Active  match  increased  primarily  because  the  total  coupling  was 
insufficient  to  cancel  the  significant  passive  reflection 
coefficient  associated  with  the  electrically  small  notch.  The 
low  antenna  gain  is  attributed  to  the  poor  passive  match,  a 
characteristic  of  electrically  small  antennas. 

6.0  CONCLUDING  REMARKS 

The  feasibility  of  using  notch  array  elements  over  a  6/1 
frequency  band  has  been  demonstrated. 
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Figure  5.  Element  Swept  Gain  (18  to  19  GHz) 


Figure  d.  Element  Swept  Gain(s)  Relative  To 
Standard  (3.0  to  6.0  GHz) 
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Figure  10.  Element  Gain,  3  to  19  GHz 
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Figure  11.  Elevation-Plane  Half-Power  Beamwidth 
(2  to  19  GHz) 


■lot 


Step  (Notch) 


Notes:  1.  Shaded  area  and  line  feed 
are  1  oz  copper  on  1/16  TK 
duroid  (r~2.2) 


2.  a  =  kb,  select  ‘K’ 

empirically  for  best  VSWR 
bandwidth 


Figure  12.  Button-Hook-Fed  Coupling  to  Slot 


276 


AN  ORBITER-ACTS  COMMUNICATIONS  EXPERIMENT 


G.  D.  Arndt,  L.  Leopold,  S.  H.  H.  Chen  -  NASA/Johnson  Space  Center 
V.  C.  Loh,  Veng  S.  Kuo,  R.  Shaw,  J.  Carl  -  LEMSCO 


Abstract 

A  proposed  communication  experiment  between  the  Shuttle  Orbiter 
and  the  ACTS  (Advanced  Communications  Technology  Satellite)  is 
described.  The  experiment  includes  20/30  GHz  technology,  a  high 
gain  antenna  (52  dB)  aboard  the  Shuttle,  and  two  ground  stations. 
Introduction 

The  Advanced  Communications  Technology  Satellite  (ACTS)  is  an 
experimental,  geosynchronous  satellite  scheduled  to  be  launched  in 
1991.  This  satellite,  transmitting  at  20  GHz  and  receiving  at  30 
GHz,  provides  spot  beams  to  fixed  ground  locations  within  the  United 
States.  It  also  has  a  program  steerable,  1  meter  antenna  which  can 
communicate  with  spacecraft  in  low-earth  orbit. 

This  paper  describes  a  proposed  communication  experiment 
between  the  Shuttle  Orbiter  and  the  ACTS.  The  objectives  of  the 
experiment  include  (1)  operate  a  high  gain  antenna  (52  dB)  with  both 
autotracking  and  computer  pointing  in  the  dynamic  environment  of  a 
maneuvering  Orbiter,  (2)  develop  and  test  space  qualified  antenna 
and  RF  technology  at  20/ 3U  GHz,  and  (3)  evaluate  multi  beam 
communication  links  for  direct  distribution  of  high  rate  data  for 
manned  space  applications.  The  experiments  will  involve  the 
Orbiter,  the  ACTS,  the  Lewis  Research  Center  Ground  Station  and  the 
Johnson  Space  Center  (JSC)  Communications  Test  Facility. 
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Initial  in-house  studies  have  been  funded  using  the  Johnson 
Space  Center's  Director's  Ui screti onary  Funds  for  fiscal  years  FV8/ 
and  FY88.  These  studies  included  system  sizing  tradeoffs.  Shuttle 
compatibility  issues,  20  GHz  electronic  fabrication  and  testing, 
equipment/technology  surveys,  and  costing.  The  study  results  are 
summarized  in  this  paper. 

ACTS  Configuration 

There  are  two  possible  modes  of  operation  for  the  ACTS  system. 
The  first  is  the  microwave  switch  matrix  mode  which  provides  a 
flexible  format  capabilty  (Figure  1).  A  ground  user  transmits  at  30 
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Figure  1.  Microwave  Switch  Matrix  Mode 


GHz  to  the  satellite  receiver,  through  the  microwave  switch,  to  the 
20  GHz  transmitters,  and  then  down  to  a  ground  station.  The  matrix 
switch  in  ACTS  is  a  4*4  switching  system  with  one  redundant  channel. 
This  makes  a  3*3  matrix  switch  available  to  the  user.  The  switch 
can  handle  up  to  three  distinct  antenna  inputs  at  any  one  time, 
regardless  of  the  signal  origin.  Three  different  signals  can  come 
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from  any  of  the  three  fixed  beams,  the  East  scanning  beam,  and/or 
the  West  scanning  beam.  Although  the  scanning  beam  can  be  in  only 
one  place  at  a  time,  it  can  move  to  different  locations  within  the 
sector.  The  beam  can  hop  around  continuously  -  up  to  4U  times  in  a 
given  frame.  Basically,  this  is  a  bent-pipe  mode  and  any  modulation 
type  in  either  a  TDMA  or  F DMA  format  can  be  supported.  The 
frequency  plan  accommodates  900  MHz  of  useable  bandwidth,  permitting 
the  transmission  of  a  very  high  burst  rate  signal,  as  well  as  a  very 
low  burst  rate  signal . 

The  second  method  of  operation  is  a  baseband  processor  mode  in 
which  the  signal  is  detected,  processed  at  baseband,  and  then 
remodulated  onto  the  downlink  carrier.  Because  the  proposed 
experiment  is  Space  Shuttle  based,  the  relative  velocity  between  the 
transmitting  and  receiving  antennas  varies  drastically,  thereby 
causing  large  doppler  shifts.  Figure  2  shows  the  doppler  profile 
for  the  proposed  experiment.  The  maximum  doppler  shift  during  the 
experiment  is  approximately  700  KHz  which  exceeds  the  operating 
capability  of  the  baseband  processor. 

Figure  2. 

Doppler  (SSO-ACTS  worst  cast  @  29.5  GHz) 
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Experiment  Configuration 


Because  of  frequency  spectrum  congestion,  future  high  capacity 
communication  links  are  moving  towards  Ka-band  (or  higher) 
frequencies.  A  high  frequency,  high  data  rate  link  to  a  low-earth- 
orbit  (LEO)  spacecraft  presents  unique  requirements  not  seen  in 
previous  programs.  This  experiment  is  intended  to  provide  answers 
to  the  following  questions: 

(1)  How  to  acquire  and  track  very  narrow  beams  of  RF  signals 
at  20/30  GHz  in  the  dynamic  environment  of  an  orbiting  vehicle?  The 
dynamic  environment  consists  of  position  changes  as  the  vehicle 
moves  along  its  trajectory  as  well  as  attitude  changes  when  reaction 
control  system  (RCS)  jets  are  fired  to  maintain  a  given  attitude. 

The  obvious  platform  to  carry  out  this  test  is  the  Space  Shuttle 
Orbiter,  which  has  an  orbital  velocity  of  about  7.b  Km/sec  and  a 
maximum  attitude  rate  of  9°/second. 

(2)  What  is  the  technology  required  to  provide  efficient  high 
gain  antennas  that  will  perform  in  this  environment?  Some  specific 
concerns  are  surface  tolerances,  feed  design,  and  tracking  mechanism 
needed  to  provide  the  tight  pointing  accuracy  required. 

(3)  How  readily  available  is  the  RF  technology  and  what  kind 
of  risk  is  involved  to  develop  space  qualified  components  such  as 
low  noise  amplifiers  at  20  GHz  and  high  power  TWTAs  at  30  GHz?  What 
is  the  noise  figure  and  DC  to  RF  efficiency  of  these  devices? 

(4)  What  kind  of  working  environment  will  the  satellite's 
multibeam  architecture  provide  for  direct  distribution  of  high  rate 
data  for  applications  such  as  the  Space  Station? 
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The  Electronics  Systems  Test  Laboratory  (ESTl)  at  JSC  will  be 
the  prime  ground  station  while  the  Lewis  Research  Center  will  serve 
as  the  secondary  ground  station.  As  described  in  an  earlier 
paragraph,  the  IF  switching  mode  of  the  ACTS  will  be  used  to  receive 
the  transmission  from  the  Orbiter  due  to  the  high  doppler  rate  (700 
KHz).  Figure  3  shows  the  experimental  configuration.  An  RF  signal 
from  either  one  of  the  ground  stations  will  be  received  by  the  ACTS 
2.2  meter  antenna  via  the  East  network.  The  signal  is  then  routed 
by  the  IF  switch  matrix  to  the  I  m  steerable  for  downlink  to  the 
Orbiter.  The  uplink  from  the  Orbiter  will  be  received  using  the  1  m 
steerable  and  routed  to  the  ESTL  via  the  East  network  and  the  3.3  u 
antenna.  At  predetermined  intervals,  a  different  antenna  beam  of 
the  ACTS  is  selected,  such  that  the  LeRC  becomes  the  ground  station 
instead  of  ESTL. 
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The  experiment  is  designed  for  the  maximum  data  rate  possible  under 
the  constraint  that  the  antenna  onboard  the  Urbiter  must  be  small 
enough  to  fit  inside  the  payload  bay.  The  system  parameters  of  the 
ACTS  are  summarized  in  the  following  tables: 


ACTS  TRANSMITTING 

TO  ORBITER 

TO  ESTL 

TO  LeRC 

Diameter 

1  m 

3.  3  m 

3.  3  m 

Gain 

such  42  dB 

50  dB 

51  dB 

TWTA  Transmitter 

10  W/43  W 

Same 

Same 

ACTS  RECEIVING 

From  Orbiter 

From  ESTL 

From  LeRC 

Diameter 

1  m 

2.2  m 

2.2  m 

Gain 

45  dB 

50  dB 

51  dB 

LNA  Noise  Figure 

5  dB 

Same 

Same 

Wich  these  parameters,  the  goal  is  to  design  a  system  that  is  able 
to  support  slow  scan  TV  (5  MbPS)  for  the  uplink  and  full  motion 
color  TV  (50  Mbps)  for  the  downlink  during  the  test.  A  more 
detailed  configuration  is  shown  in  Figure  4  where  ESTL  will  be  the 
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ORBITER 


Figure  4.  Detail  Experiment  Setup 

experiment  operations  as  well  as  the  data  analysis  center.  Commands 
and  telemetry  are  routed  to  the  Grbiter  via  the  S-band  or  Ku-band 
system  and  coordination  with  LeRC  will  be  by  prior  arrangement  and 
real  time  audio  circuits.  System  configurations  for  each  of  the 
ground  stations  are  shown  in  Figure  5.  The  uplink  data  from  ESTL 
will  be  either  5  Mbps  of  PN  code  or  digitized  slow  scan  TV  which 
will  be  BPSK  modulated  onto  the  29.6  GHz  carrier.  The  ESTL  power 
amplifier  will  be  a  1U  W  TWTA  and  the  antenna  will  be  a  1.8  m 
reflector  with  monopulse  feed.  The  ESTL  ground  receiver  will 
receive  a  rate  1/2,  convol uti ona 1 ly  encoded  55  Mops  (110  Mbps  total) 
signal  at  19.7  GHz.  This  unbalanced  QPSK  signal  will  be  demodulated 
into  the  1/ Q  components .  If  PN  spauences  are  transmitted,  bit  error 
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rate  measurements  will  be  performed  while  the  uplink  and  downlink  RF 
powers  are  varied.  If  a  TV  signal  is  transmitted,  the  base  band 
data  will  be  routed  to  the  TV  laboratory  at  the  Tracking  and 
Communication  Division  at  JSC  for  evaluation  of  future  EVA/Space 
Station  applications.  The  LeRC  station  will  be  an  adjunct  to  the 
ESTL  site  and  will  mainly  be  used  to  evaluate  the  effect  of  real 
time  routing  of  the  110  Mbps  data  by  ACTS  switching  network.  Its 
equipment  will  consist  of  an  uplink  modem  using  a  5  Mbps  PN 
generator  as  the  data  source  while  the  receiver  will  consist  of  a 
QPSK  demodulator  and  a  high  data  rate  recorder  where  the  data  will 
be  saved  for  post  analysis.  The  BPSK  signal  sent  from  either  of  the 
ground  stations  will  be  demodulated  in  the  Orbiter  receiver  and  then 
turned  around  to  modulate  the  Q  channel  of  the  downlink  QPSK  signal. 
If  slow  scan  TV  is  used  as  the  uplink,  the  signal  will  be  converted 
to  the  Orbiter  CCTV  format  and  distributed  via  the  CCTV  network. 

The  BO  Mbps  high  rate  data  will  be  generated  either  from  a  PN 
generator  or  via  a  special  TV  signal  digitizer.  The  TV  signal 
obtained  from  the  Urbiter  CCTV  system  allows  the  crew  to  select  any 
of  the  cameras  located  in  the  payload  bay.  This  signal  will  then  be 
used  to  modulate  one  channel  of  the  modulator. 

The  antenna  used  in  the  Orbiter  is  the  same  1.8  m  antenna  used 
in  ESTL  except  that  it  will  be  mounted  inside  the  payload  bay  for 
safety  reasons  as  shown  in  Figure  6.  Although  the  scan  angle  is 
limited  and  the  pattern  is  affected  by  the  bulkheads  of  the  payload 
bay,  the  cost  is  minimized  by  not  mounting  the  antenna  on  a  moveable 
boom. 
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Figure  5:  Ground  Strut  ions  and  Orb  iter  Systems 


With  the  links  configured  as  shown  in  Figures  4  through  6,  the 
performance  is  satisfctory  between  all  the  stations  and  the  Orbiter. 
Using  a  50  watt  TWTA,  30  GHz  transmitter  on  the  Orbiter,  the  worst 
case  downlink  channel  (50  Mbps)  has  5.7  dB  of  RF  circuit  margin  as 
given  in  Tables  A  and  B.  A  summary  of  the  requirements  for  the 
different  types  of  antennas  is  given  below.  Development  efforts 
will  concentrate  on  the  1.8  m  antenna,  low  noise  amplifiers,  TWTAs, 
and  the  necessary  baseband  equipment. 


SSO  to  ACTS 


ACTS  to  SSO 


Pt  (W) 

50 

43 

TX  Ant  (m) 

18 

1 

Polarization 

Circular 

Horizontal 

Pointing  Accuracy 

±0.06' 

±0  15' 

RX  Ant  (m) 

1 

1.8 

Polanzation 

Vertical 

Circular 

Pointing  Accuracy 

±0.1' 

±0  09° 

Circuit  Margin  (dB) 

57 

80 

ESTL  to  ACTS 

ACTS  to  ESTL 

Pt(W) 

10 

43 

TX  Ant  (m) 

18 

3.3 

Polanzation 

Vertical 

Horizontal 

Pointing  Accuracy 

±0.06“ 

±0.05' 

RX  Ant  (m) 

2.2 

18 

Polarization 

Vertical 

Horizontal 

Pointing  Accuracy 

±0  075“ 

±0  09“ 

Circuit  Margin  (dB) 

11.7 

120 

Hardware  Studies 

As  part  of  the  initial  definition  of  the  experiment,  several 
in-house  studies  into  hardware  availability  are  underway.  Results 
to  date  are  summarized: 
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TABLE  A 


CIRCUIT  MARGIN -SSO  to  ACTS  (50  mbps  Channel) 


Pt  (3>  50  W 
Ckt  Loss 

Polarization  Loss  (circular  to  vertical) 
Pointing  Loss  (±0.06°  SSO,  ±0.1°  ACTS) 
TX  antenna,  1.8m@  50%  r| 

EIRP 

Path  loss  (30  GHz,  22000  nmi) 

RX  antenna,  1  m  @  30%  (per  RCA) 

Ckt  Loss 
Total  Rec.  Power 

QPSK  modulation  loss  (80  %) 

QPSK  demodulation  loss 

Ts  (5  dB  NF,  Ta  =  290°  K) 

No 

Bit  rate  bandwidth  (50  Mbps) 

^b^o  in  bit  rate  bandwidth 

Req'd  Eb/NQ  (10'5  BER) 

Coding  gain  (r=1/2,  k=7) 

Bit  sync  Loss 
Actual  Eb/NQ  required 


17.0 

dBW 

-2.0 

dB 

-3.0 

dB 

-0.5 

dB 

52.0 

dB 

63.5 

dBW 

-214.2 

dB 

45.0 

dB 

-3.0 

dB 

-108.7 

dBW 

-1.0 

dB 

-1.0 

dB 

29.6 

dBK 

-199.0 

dB(W/Hz) 

77.0 

dBHz 

11.3 

dB 

9.6 

dB 

5.0 

dB 

-1.0 

dB 

5.6 

dB 

5.7 

dB 

Margin 


TABLE  B 


CIRCUIT  ...MARGIN  -  ACTS  to  ESTL  f50  mbps  Channph 


Pt  @  43  W 

16.3 

dBW 

Ckt  Loss 

-3.0 

dB 

Polarization  Loss  (both  are  horizontal) 

-1.0 

dB 

Pointing  Loss  (±0.05°  ACTS,  ±0.09°  ESTL) 

-0.5 

dB 

TX  antenna,  3.3  m  fix  mounted  on  ACTS 

50.0 

dB 

EIRP 

61.8 

dBW 

Path  loss  (20  GHz,  22000  nmi) 

-210.7 

dB 

RX  antenna,  1 ,8m(5)  50%  r|  at  ESTL 

48.5 

dB 

Ckt  Loss 

-2.0 

dB 

Total  Rec.  Power 

-102.4 

dBW 

QPSK  modulation  loss  (80  %) 

-1.0 

dB 

QPSK  demodulation  loss 

-1.0 

dB 

Ts  (5  dB  NF,  Ta  =  290°  K) 

29.6 

dBK 

No 

-199.0 

dB(W/Hz) 

Bit  rate  bandwidth  (50  Mbps) 

77.0 

dBHz 

Eb/NQ  in  bit  rate  bandwidth 

17.6 

dB 

Req'd  Eb/NQ(10*5  BER) 

9.6 

dB 

Coding  gain  (r=1/2,  k=7) 

5.0 

dB 

Bit  sync  Loss 

-1.0 

dB 

Actual  Eb/NQ  required 

5.6 

dB 

Margin 

12.0 

dB 
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(1)  Low-Noise  Amplifier  -  One  of  the  key  components  in  the 
implementation  of  the  ACTS  experiment  is  the  receiver  low-noise 
amplifier.  This  low-noise  amplifier  should  add  minimum  noise  power 
to  the  received  K-band  signal  for  the  system  to  achieve  maximum 
signal -to-noise  ratio.  In  addition,  the  gain  of  the  LNA  should  be 
high  enough  to  mask  second  stage  noise  contributions  of  the  filters, 
mixers,  and  pre-ampl i fiers . 

Due  to  the  recent  advances  in  HEMT  (high  electron  mobility 
transitor)  device  technology,  development  of  a  prototype  hybrid 
microwave  integrated  amplifier  has  been  undertaken  at  NASA/ JSC. 

This  amplifier  has  been  implemented  using  a  0.5  micron  gate  length 
HEMT  device.  The  design  goals  for  the  single-stage  amplifier  unit 
are  a  gain  of  7  dB  and  a  noise  figure  of  4.5  dB.  The  amplifiers 
have  been  realized  in  a  single-ended  confi guration  with  maximum 
power  gain  matching  networks  on  the  input  and  output.  The  complete 
low-noise  amplifier  is  a  3-stage  cascaded  circuit  with  interstage 
gain  compensation  networks  to  achieve  20  +/-  1  dB  gain  and  5.0  dB 
noise  figure. 

(2)  Pan-Tilt  Units  vs.  Gimbals  -  Gimbals  are  normally  used  for 
closed  loop  tracking  and  positioning  applications  of  high  gain, 
spaceflight  antennas.  After  being  space  qualified,  gimbals  are  very 
reliable  and  efficient  during  acquisition  and  tracking  manuevers. 
However,  the  cost  can  reach  $3M  to  $8M  for  design,  development, 
testing,  and  space  qualifying  of  the  gimbals  and  associated  control 
el ectronics . 
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A  possible  alternative  to  the  gimbals  are  pan-tilt  devices 
employing  stepper  motors.  The  pan-tilt  unit  makes  use  of  a  rotary 
incremental  activator  to  produce  a  simple  and  flexible  way  of 
gimballing  a  load.  The  rotary  incremental  activator  is  a  compact, 
closely  integrated  design  made  up  of  two  key  elements  -  a  motor  and 
a  harmonic  drive.  These  activators  have  been  used  in  such  varied 
applications  as  antenna  deployment  mechanisms,  latching  devices, 
large  mirror  positioning  systems,  antenna  gimbal  drives,  and 
instrument  pointing. 

If  a  pan-tilt  unit  can  operate  efficiently  for  monopulse 
tracking  applications  of  large  space  antennas,  the  cost  advantage 
over  a  conventional  gimbal  system  makes  it  an  attractive  alternative 
for  this  experiment.  A  detailed  engineering  comparison  of  the  two 
techniques  is  underway. 

Urbiter  Antenna  Design  Considerations 

The  requirement  for  a  large  Orbiter  antenna  at  20/30  GHz  may 
conflict  with  the  general  Orbiter  requirement  for  a  light  weight 
design.  A  simulation  program  was  developed  to  ascertain  the 
requirements  for  rigidity,  smoothness,  and  mechanical  precision  in 
the  1.8  meter  reflector  and  feed  in  order  that  50  dB  could  be 
reliably  obtained  from  the  antenna.  It  is  assumed  that  the  antenna 
uses  a  four  horn  monopulse  autotrack  feed  with  the  feed  cluster 
located  directly  above  the  reflector's  vertex.  For  the  sake  of 
compactness,  a  shallow  F/D  of  U. 35  was  assumed.  The  aperture 
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9 

illumination  AI  (il)  is  given  by:  AI  (a)  =  1  -  ("v^m)  where  /t 
is  the  local  radius  and^m  is  the  maximum  radius.  An  overall 
efficiency  of  40%  is  assumed. 


Table  C 

Effects  of  Reflector  Deformation  on  Antenna  Performance 


Deformation 

Type 

D(mm) 

Effect  on  Radiation  Pattern 

Type  1 

2 

Loss 

Beam 

of  Boresight  bain  =  0.22  dB 
Shift  c  o.U7°  (in  AZ) 

4 

Loss 

Beam 

of  Boresiynt  Gain  =  0.95  dB 
Shift  *  0.14°  (in  AZ) 

1  y  pe  2 

/7 

2 

Loss 

Beam 

of  Boresight  Gain  =  0.19  dB 
Shift  =  0.0° 

\\ 

4 

Loss 

beam 

of  Boresight  Gain  =  1.34  dB 
Shift  *  0.0° 

Type  3 

2 

Loss 

Beam 

of  Boresight  Gain  =  1.68  dB 
Shift  X  0.2°  (in  AZ ) 

C  - 

4 

Loss 

Beam 

of  Boresight  Gain  =  7.9  dB 
Shift  £  0.4°  (in  AZ) 

Type  4 

D- 

2 

Loss 

Beam 

of  Boresight  Gain  =  1.36  dB 

Shift  =  0.0 

K 

4 

™  Loss 
Beam 

of  Boresight  Gain  =  7.39  dB 
Shift  -  0.0 
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a.  Reflector  Deformations 

Reflector  deformations  may  be  caused  by  thermal, 
mechanical,  or  inertial  stresses,  or  by  distortion  from  a  true 
parabola  during  the  fabrication  process.  Representative 
deformations  and  the  associated  degradations  are  summarized  in 
Table  C.  A  typical  gain  versus  degradation  plot  is  shown  in  Figure 
7. 

Types  2  and  4  are  radially  symmetrical  deformations  and 
consequently,  produce  gain  loss  but  no  beam  shift.  Types  1  and  3 
produce  both  beam  shift  and  gain  loss.  The  results  indicate  a 
reasonable  distortion  limit  of  2  mm  produces  .2  to  1.7  db  of 
antenna  loss,  depending  upon  the  particular  type  of  deformation. 

b.  Surface  Roughness 

Figures  8  and  9  indicate  that  the  surface  of  the  reflector 
should  be  smooth  to  within  approximately  one-half  millimeter  (1mm) 
in  order  to  keep  the  gain  degradation  to  0.3  db  and  the  first  side 
lobe  level  increase  to  1.5  dB.  The  surface  material  and 
fabrication  technique  must,  therefore,  be  capable  of  obtaining  and 
maintaining  smoothness  to  0.5  mm  or  less. 

c.  Lateral  Feed  Mislocation 

Laterial  feed  mislocation  is  defined  here  as  a  mislocation 
along  any  line  orthogonal  to  the  boresight  axis  of  the  reflector. 
Displacement  is  in  the  "+y"  or  "+Azimuth"  direction.  For  small 
lateral  feed  mislocation,  the  principle  effect  is  beam  steering  as 
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GAIN  VS.  SURFACE  ROUGHNESS 


shown  in  Figure  1U  with  little  loss  of  peak  gain.  However,  loss  of 
boresight  gain  can  result  as  can  be  seen  from  Figure  11. 
d.  Boresight  Axis  Feed  Mislocation 

Feed  mislocation  "in"  or  "out"  along  the  reflector 
boresight  axis  causes  defocusing  of  the  antenna.  This  results  in  a 
loss  of  gain  and  a  smearing  of  the  first  sidelobe  into  the  main 
beam.  (It  is  interesting  to  note  that  the  smearing  effect  starts 
to  occur  before  serious  gain  loss.  This  phenomenon  could  be  used 
to  avoid  false  tracking  on  the  first  side  lobe.)  This  type  of  feed 
mislocation  does  not  cause  the  beam  to  shift  from  boresight. 

Figure  12  shows  the  relationships  between  boresight  axis  feed 
mislocation  and  gain. 

Summary; 

A  20/30  GHz  experiment  between  the  ACTS  and  the  Orbiter  is 
technically  feasible  with  much  of  the  spacecraft  hardware  derived 
from  ground  equipment  already  developed  as  part  of  the  ACTS 
program.  The  experiment  can  be  self-contained,  with  no  scarring  to 
the  Orbiter.  In  addition  to  the  space  flight  aspect  of  operating  a 
high  gain,  high  rate  data  system  from  the  Orbiter,  the  advantages 
of  direct  satellite  broadcast  to  the  Johnson  Space  Center  can  be 
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MOM  Analysis  of  a  Finite  Length  Slot 
In  an  Infinite  Conducting  Sheet 


Henry  A.  Karwacki 

Sanders  Associates,  Inc. 
Nashua,  New  Hampshire,  03061 


Abstract 


This  paper  presents  a  method  of  moments  analysis 
of  a  finite  length  and  width  slot  in  an  infinite 
conducting  sheet. 

The  analysis  uses  entire  domain  basis  and  testing 
functions  (i.e.  sinusoids)  and  image  theory  to  model 
the  scattering  from  a  finite  length  and  width  slot  in 
a  ground  plane  a  for  an  incident  TM  wave. 

The  computations  performed  for  a  single  cosine 
basis  function  and  for  a  five  term  summation 
demonstrated  very  little  difference  since  all  higher 
order  terms  were  9  dB  or  more  below  the  first  term. 

The  solutions  also  demonstrated  that  the  form  of  the 
scattered  far  field  was  the  same  as  the  far  field  of  a 
radiating  rectangular  aperture. 
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Introduction 


1  . 

The  scattered  field  from  a  narrow  width,  finite 
length  slot  in  an  infinite  perfectly  conducting  sheet 
excited  by  an  incident  transverse  magnetic  (TM)  plane 
wave  (as  shown  in  Figure  1)  is  examined.  A  method  of 
moments  solution  is  formulated  using  entire  domain 
basis  functions  with  a  Galerkin  testing  procedure.  It 
is  of  interest  to  quantitatively  compare  the  relative 
amplitudes  of  the  higher  order  slot  mode  excitations 
to  the  lowest  order  slot  mode.  It  is  evident  that  the 
entire  domain  expansion  needed  converges  very  rapidly. 

The  inner  products  which  fill  the  system  matrix 
usually  involve  two  sets  of  integrations,  in  this 
formulation  the  inner  products  are  reduced  to  a  single 
set  of  numerical  integrations  thereby  conserving 
computer  resources.  If  the  field  produced  by  the  slot 
magnetic  current  is  considered  as  the  convolution  of 
the  current  with  the  free  space  Green's  function,  it 
is  possible  to  interchange  the  testing  function  and 
the  Green  's  operator  through  a  proper  change  of 
variables.  The  Green's  operator  will  then  act  on  the 
convolution  of  the  currents  and  the  testing  functions. 
This  formulation  can  also  be  used  to  compute  the  self- 
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impedance  of  a  dipole  or  slot  or  the  mutual  impedance 
between  two  dipoles^ 


2.  Formulation 

Consider  a  slot  which  is  oriented  as  shown  in 
Figure  1  on  a  perfectly  conducting  sheet  of  infinite 
extent  in  the  x  and  y  dimensions  and  vanishingly  thin 
in  the  z  dimension.  The  slot  has  a  finite  length  L  and 
a  finite  width  W.  A  plane  wave  which  is  TM  to  the  slot 
is  incident  on  the  sheet.  The  incident  H-field  is 
given  by 


H>  = 
me 


H . 
me 


(1  ) 


The  boundary  conditions  which  must  be  met  at  the 
slot/sheet  interface  is  that  the  tangential  E-field 
and  the  total  H-field  across  the  slot  is  continuous, 
and  the  tangential  E-field  on  the  sheet  is  zero.  The 
tangential  electric  field  in  the  slot  excited  by  the 
plane  wave  can  be  replaced  by  equivalent  magnetic 
current  sources  distributed  over  the  slot  area  on  both 
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sides  of  the  conductor  after  closing  the  slot  with  a 

conductor.  By  applying  image  theory  as  shown  in  Figure 

2,  the  field  on  the  incident  side  of  the  sheet  is 

reduced  to  the  field  produced  by  twice  the  magnetic 

current  in  free  space,  the  incident  field  and  its 
r  3  4  ] 

image  '  .  Setting  the  total  magnetic  field  at  the 

slot  on  side  1  equal  to  that  on  side  2  an  equation  is 
obtained  relating  the  incident  field  to  the  magnetic 
field  produced  by  the  magnetic  current  source  M  only. 


H(M)  =  H 


inc 


(2) 


The  electric  vector  potential  derived  from  the 
slot  current  is  given  by 


F>  = 


47 r 


M 


w/2  ,  L/2 


-w/2  -L/2 


- jk | r-r  ' 
dr' 


r-r 


M>=Mx 


/  2  2  2 
- jk\/(x-x  ' )  +  (y-y  ' )  +  ( z-z  ') 

M  x  — y - 2 - 2 - ~~2 

\/( (x-x  ' )  +  (y-y  ' )  +  (  z-z  ' ) 

(  3 


and  the  electric  and  magnetic  fields  are  given  by 


1 

j^e0 


V>x  H> 


F> 


(5) 


(see  appendix  A  for  definitions  of  variables) 

The  basis  functions  chosen  for  the  method  of 
moments  solution  are  entire  domain  basis  functions 
based  upon  a  sum  of  cosine  functions  of  integral 
multiples  of  x  having  support  over  the  slot  region. 

The  testing  functions  are  chosen  to  be  the  same  as  the 
basis  functions  according  to  Galerkin's  method. 

The  equivalent  magnetic  current  in  the  slot  can 
be  expressed  as 


M> ( x  ' , y  ' , z ' )  = 

-  N 

x  )  a  cos  (  ^7—  x  ' )  •  [  u  ( x  '  +  L/2  ) -u  ( x  '-L/2  )  ] 
n=l 

•  [  u  ( y  '+w/2  ) -u  (  y  '-w/2  )  }«  ^  ( z ') 

(6) 

which  when  substituted  into  the  equation  for  the 
electric  vector  potential  F  yields 


,  w/2^  L/2 


F>=  x 


47T 


L 


N 


—  00  -w/2  —L/2 


n  =  l 


-jkR 

a  cos(~x');  $  (  z  '  )dx  'dy  'dz 
n  l  h 

(7 


30] 


where ; 


r=V'  x_x  ' > 2+<y-y  ' ) 2+ ( z-z  ' ) 2 
then 


>-oi  r— 2n  ^  - jkR 

=  x  k‘^1  }_  ancos(^x')|  S(z')dx'dy'dz 

J  s .. J  T  ,^n=l 


-oo  -w/2  -L/2 


nw/2  ✓  L/2 

J  V2[L  anCOS  (  '3l  s  (  ^  )dx  'dy  'dz  ' 

_  __/o  ^  n  x 


-OO  -w/2  -L/2 


This  equation  can  be  rewritten  as  convolution 


integral  in  x 


P(x)  =  cos(n^~ )  [u(x+  -  u(x-  ]  ; 


Q  (  x )  = 


-jkR 


then 


F  (  x  )  ©  Q  (  x  )  dy  '  £  (  z  '  )  dz' 


-w/2 


+>  a. 


nw/2 

V>(V'(P(x)  ©  Q(x)  )dy  £(z  ')dz  ' 
_ oo  -w/2  i 


( where  ©  represents  convolution) 
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Applying  the  distributive  properties  of  the 
convolution  integral  with  regard  to  differentiation, 
the  Laplacian  can  operate  on  either  P  or  Q  prior  to 
performing  the  convolution. 


V2(P(x)0Q(x)  =  P(x)  0  (V  2Q(x)  )  =  (  V2P(x)  )  0  Q(x) 

Since  the  incident  magnetic  field  has  only  an  x 
component  the  Laplacian  reduces  to  the  second  partial 
derivative  on  x. 


2  2 

2  P(X)  =  -(—£-)  cos(-^-x)  [u(x+  -r- )  -  u(x-  <-y)  ] 
dx  u  z 

+  ^^[sin(^)[S(x+  —)  -S(x-  -—)]]  (10) 

+  cos(£y)  [  g'(xv  -±-)  -  $'(x-  -i-)] 

which  for  n  even  reduces  to; 

2 

-  (~y)  cos(~)[u(x+  — J- )  -  u(x-  **2~)] 

+  cos(~-)  [g'(x+  -~)  -  S'(X-  4p)] 
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and  for  n  odd  ; 


2  2 

---  P  ( X )  =  COS  (  ““X  )  [  u  ( x+  —7- )  -  u(x-  -4-)] 

C  ij  L  Z  Z 

dx 


+  -7-^  sin(^yL )  [  $(x4  )  +  S(x-  )  ] 

L  Z  2  *  (12) 


The  choice  of  n  even  implies  that  an  E-field 
maximum  at  x  =  +/-  L/2  which  is  inconsistent  with  the 
physical  assumptions,  therefore  the  choice  of  n  odd  is 
employed.  Substituting  (11)  into  (8)  results  in  the 
following  equation  for  the  magnetic  field  in  the  x 
direction 


N 


n 


r  w/2  ,  L/2 
-w/2  -L/2 


L  )  )cos(“-x  ') 


-  jkR  ' 

e  J  . 

T  dx  d* 


n  r w/2 

-LV-  J 

n=1  -w/2 


(  2n7f  2L  '  ,  n'TT'  e _ 

(  t"  -  ~5sir‘(  — ^  dV 


n  If 


( 1  J) 


where  R'  is  R  evaluated  at  z  =  0  , 


R 


=  \Zx-x 


>  \  2  ,  ,.2 

)  +  (y-y  ) 
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and  R''  is  R  evaluated  at  z  =  0,  x'=  L/2 


R"  =  \/(x*  L/2  )  2+  (  y-y  ' )  2 


The  choice  of  testing  function  is  made  according 
to  Galerkin's  method^  Hence  an  element  of  the  set 
is  given  by 


T  (  x ) 

m  odd 


cos(^x)  [u(x+  -  u(x-  -^-)]  • 


{ u  ( y  +  -tt)  -  u(y-  -~)]  (14) 

When  (2)  is  tested  by  (13)  over  the  support  of 

the  slot  a  set  of  simultaneous  linear  equations  result 

from  which  the  unknown  current  coefficients  a  can  be 

n 

determined . 


<Tj (x)  ,  H^(x)>  -  - 


-  -  <Tj(x)  ,  Hn(x)>  aj 


<Tn(x)  ,  H^(x)> 


-  -  <Tn(x)  ,  HN(x)>|  aN 


<Tl(x)  ,  Hj(x)> 


<Tn(x)  f  Hn(x)> 


an  arbitrary  inner  product  is  evaluated  as 


<Tm(x>  , 


HS>  - 
X 


L/2  , w/2 


2 

(k  -( 


-L/2  -w/2 


,  n  <7T  -  >  , 

cos(~x  )cos( 


-  jkR  ' 

— dx  'dy  'dx 
K 


L/2  /•  w/2 


—  -i  ^ 

-|P  -  ^)sin(Brlcos(!!ax)S-,.dy'dx 


-L/2  -w/2 


L/2/'  w/2 


-L/2  -w/2 


,2x1^  2L  \  -n'T,  ,m'7T  \G 

(16) 


The  above  equations  can  be  reduced  and  simplified 
by  further  application  of  convolution  methods  and 
manipulation  of  variables  to  yield 
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4Lm 


<T  (x)  ,  H  >  = 
m  x 


ENa 

n=l  n 


L  /  w/2 

( 

0  -w/2 


(m2-n2)fl~ 


[ sin ( ( L.x  - )  ,cos ,  iffl±ai2Cx  - , 


+  sin  (  ( L- x  - ) )  cos  ( - ) ) .  < k2- ( sac ,  I 


e-jk\/  x'2  +  y'2 

*  ~r~r  .2  dii'dx' 

\/  X  +  y 


L/2  w/2 


-L/2  -w/2 


2n 


‘  ^)sin(!1f)cos(Sfx: 


,  m  '77'  ,  e 

COS  (  — “'X  )  — 


V  (x-1  ""2‘ - '  ‘  2 


■1/2  j  +  (  y-y 


/  T  T  L/2  w/2 

-jk\/(x+L/2)^+(y-y V 
“7 - 5 - 2  dy  dx  ■* 

\/ (x+l/2)/+(y-y')z 

-L/2  -w/2 

- jk\/(x-L/2 ) 2+ (y-y  ' ) 2 


/  / 


2n?T  2L^  j  •  ,  2l2C_. 

L  n^T)sin(  2  ’ 


dy  'dx 


(17) 


and 


<T  (x)  ,  H1>  = 

m  x 


{  L/2 
-L/2 


,rr V7C  .  -ikx  sin©.cos$.  , 
cos  (  — j—  x)e  -  1  1  dx 


(18) 
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[-2cos(  K  ) -^-sin{ +  2k  sin  (  If  )  sin©.  cos$  .  cos  ( )  ] 
—  -Li  L  _  O  _  1  1  Z 

,,  2  .  2Q  2j  .  ,  ir,  -Tf  .  2 

(k  sm  8. cos  0.)  -  {--L) 
oil  L 

where  =  sine^cos^ 

4.  Computations 

The  solution  of  the  matrix  equation  was  performed 
using  numerical  integration  and  matrix  inversion  for  N 
=  1  ,  3  ,  5  ,  7  ,  and  9.  The  far  electric  field  was 
computed  for  each  "mode"  and  the  sum  of  these  first 
five  odd  modes.  The  results  are  plotted  for  phi  =  0° 
and  90°  as  a  function  of  the  far  field  angle  theta  (©) 
in  Figures  3  through  12  for  each  of  the  modes  and 
Figures  13  and  14  for  the  sum  of  the  five  modes.  The 
plots  of  Figures  15  and  16  are  for  phi  =  0°  and  90°  as 
a  function  of  the  angle  of  incidence  of  the  plane  wave 
excitation.  It  is  significant  to  note  the  following; 
the  computed  far  field  for  the  sum  of  the  first  five 
odd  modes  and  that  for  the  first  mode  N  =  1  are  almost 
identical,  the  relative  amplitude  of  the  higher  order 
mode  fields  are  at  least  9  Db  or  more  below  the  first 
order  mode  field,  and  the  form  of  the  scattered  far 
field  is  the  same  as  would  be  expected  from  a 
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rectangular  aperture.  It  is  also  of  note  that  for  a 
narrow  slot  (w<<^.)  there  can  be  no  transverse 
electric  wave. 

5.  Conclusions  and  Future  Application 

The  method  yields  the  scattered  field  from  a 
transverse  magnetic  wave  incident  on  finite  length  and 
width  slot  in  an  infinite  ground  plane  even  for  a 
relatively  small  number  of  terms  but  the  computations 
must  be  validated  by  comparison  to  measured  data.  The 
future  application  of  this  formulation  will  involve 
the  solution  for  a  transverse  electric  incident  wave, 
superposition  of  the  TM  and  TE  solutions  and  the 
extension  to  a  solution  for  a  pair  of  orthogonal 
slots . 
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APPENDIX  A 


Variable  Definitions 

a  n  th  coefficient  of  the  basis  function  sum 
from  1  to  N 

F>  electric  vector  potential 

H>,E>  magnetic  and  electric  field  vectors 

2 

k  wave  number  

L  slot  length 

-> 

M  magnetic  current  vector 

n  summation  index 

P(x),  Q(x-x')  functions  of  x,  x-x' 

R,  r  radius 

T(x)  testing  function 

w  width  of  slot 

theta  incident  polar  angle 
phi  incident  polar  angle 
inc,  i  subscript  for  incident  quantity 
scat,  s  subscript  for  scattered  quantity 

u[x'+_  -?j-]  unit  step  function  1  for  x'  >  *  -|- 
0  elsewhere 

(z')  impulse  function  1  when  z'  =  0  ,  0  elsewhere 
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NEAR-FIELD  BISTATIC  RCS  rEASURETlENT  AT  BDM 


R.  Rogers  and  E.  Farr 
The  BDM  Corpora t ion 
1801  Randolph  Road  SE 
Albuquerque,  NM  87106 


ABSTRACT 

The  techniques  of  near-field  antenna  pattern  measurement  can 
be  extended  to  near-field  RCS  measurerrent.  The  motivation  for 
doing  so  is  precisely  the  same  as  that  for  near-field  antenna 
measurements ;  i.e.,  the  convenience  of  an  indoor  antenna  range, 
and  an  improvement  in  accuracy.  Although  the  near-field 
neasurenent  problem  is  solvable  in  principle  in  a  manner  analogous 
to  the  neai — field  antenna  problem,  it  requires  a  significantly 
larger  amount  of  tine  to  take  the  necessary  data,  and  to 
subsequently  process  the  data  to  obtain  useful  quantities. 

BDM  is  currently  involved  in  an  on-going  program  to  evaluate 
the  feasibility  of  near-field  bistatic  RCS  measurements.  At  the 
time  of  this  viriting,  a  complete  set  of  metherratics  has  been 
formulated  to  handle  the  probe  correction  and  data  processing. 

The  hardware  has  been  built,  software  development  is  near 
completion,  and  the  analysis  of  canonical  scattering  objects  has 
been  completed.  Experimental  data  soon  to  be  taken  for  these 
objects  v/ill  be  presented.  It  is  hoped  that  the  technique  will 
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prove  to  be  a  practical  approach  to  RCS  measurements. 

1.0  INTRODUCTION 

In  this  paper  we  discuss  am  on-going  project  to  evaluate  a 
plane-wave  scattering  iratrix  (PUSH)  approach  to  measuring  bistatic 
radar  cross-section  (BRCS).  Arrays  of  measurements  obtained  with 
separate  transmit  and  receive  probe  antennas  are  computer 
processed  to  yield  the  complex  bistatic  scattering  matrix  as  a 
function  of  illumination  and  observation  angles.  The  data 
acquisition  and  computer  processing  are  an  exact  sampled  data 
implementation  of  the  continuous -domain  PWSM  theory,  provided  that 
attention  is  paid  to  bandlimit ing  and  sampling  constraints. 

A  particularly  appealing  aspect  of  this  approach  to  BRCS  is 
that  bistatic  measurements  can  be  made  at  angles  very  close  to  the 
monostatic.  A  significant  drawback  is  the  large  nunfcer  of 
measurements  that  must  be  made. 

Ue  will  discuss  in  turn  the  conti nuous -domain  PIBM  theory  of 
BRCS,  the  sampled-data  formulation  of  that  theory,  and  the 
experimental  work  in  progress. 

2.0  FLAKE-WAVE  SCATTERING  MATRIX  BRCS  THEORY 

This  work  is  based  upon  a  three-antenna  plane-wave  scattering 
matrix  theory  derived  by  Dinallc^.  The  goal  is  to  calculate  for 
an  arbitrary  object  a  (mono -frequency)  scattering  nmtnx 


322 


lies'11’1) 

Ine«(k’1) 

u«e'^’^) 

Iiw'1'’1i 

(1) 


1  =  1  e  +  1  e  +  le 

k  x  y  y  z  z 

A  a  A 

k  =  k  e  +  ke  +  k  e  . 

xx  y  y  z  z 


which  describes  the  scattering  of  a  unit  incident  plane  wave  by  an 
abject.  Incident  and  scattered  propagation  vectors  are 

A  A  A 

(2a) 
(2b) 

The  elements  of  the  scattering  rratrix  I  (k,i)  are  specified  by 
means  of  two  polarization-related  indices  and  two  propagation 
tectors.  For  exanple,  I^^(k,l)  refers  to  the  0-conponent  of  the 
angular  spectrum  of  the  wave  that  is  scattered  viien  the  object  is 
illuminated  by  the  ^-conponent  of  the  angular  spectrum  of  the 
incident  wave.  Taking  all  four  elements  together,  I  (k,l) 
specifies  how  the  9  and  components  of  the  angular  spectrum  of 
the  incident  wane  1  are  scattered  into  S  and  £  conpanents  of  the 
angular  spectrum  of  a  scattered  wave  k. 

The  scattering  of  an  arbitrary  wavefront  by  the  object 
foilavre  directly,  since  the  superposition  principle  allows  us  to 
create  complex  wavefronts  as  a  superposition  of  plane  waves. 


3.0  CALCULATING  SCATTCRING  MATRIX 

In  order  to  calculate  the  scattering  matrix,  we  will  first 
calculate  the  9-  and  tf-conponents  of  the  scattered  angular 
spectrum  over  the  area  scanned  by  the  receive  probe,  given 
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illumination  of  the  test  object  by  the  transmit  probe  in  a 
particular  orientation  and  location.  This  process  is  repeated  tor 
a  set  of  transmit  probe  locations. 

The  transmit  probe  is  then  rotated  to  illuminate  the  test 
object  with  a  different  angular  spectrum,  and  we  repeat  the 
measurenents  and  calculations  to  determine  the  8-  and  ^-conponents 
of  the  scattered  angular  spectrum  for  this  second  transmit  probe 
antenna  orientation  and  set  of  locations. 

The  laboratory  measurements  for  bistatic  near-field 
scattering  consist  of  gain  and  phase  measurements  made  by  a 
receive  probe  as  it  is  swept  through  a  pattern  of  probe  locations. 
The  receive  probe  scan  pattern  is  swept  repeatedly  (and  gain  and 
phase  measurements  made)  for  a  set  of  transmit  probe  locations,  as 
the  transmit  probe  itself  steps  through  a  set  of  locations  in  the 
transmit  probe  scan  pattern. 

The  transmit  probe  and  receive  probe  we  use  are  identical, 
although  they  need  iiot  be  so.  The  receive  probe’s  receiving 
character ist ic 


V>  - 


Wk> 


Wfe> 


(3) 


and  the  transmit  probe's  transmitting  characteristic 
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W‘> 


(4) 


2 

are  assumed  to  be  known  (either  analytically  or  enpirically  ).  As 
in  the  case  of  near-f ield  antenna  neasurements ,  neither  probe 
antenna  should  have  nulls  in  a  direction  for  which  one  wishes  to 
calculate  the  scattering  matrix. 

The  transmit  characteristic  is  related  to  the  antenna’s 
fai — field  E(r)  a s^ 


=  - Y -  limit  ^  E(t')  (5) 

2nk  cos#  r  -»  o>  '  ' 

For  a  reciprocal  antenna,  the  transmit  and  receive  characteristics 
are  related  by 


Y  cos# 

- %<-*>  !6> 

O 

where  Yq  is  the  characteristic  admittance  of  the  waveguide  feed, 
and  j7q  is  the  characteristic  admittance  of  free  space. 

The  scattering  object  is  at  the  origin  of  a  coordinate  system 
S  (Figure  1).  l*te  define  1  as  the  propagation  vector  of  the 
incident  plane  wave  (from  the  transmit  probe),  r^  as  the 
coordinates  of  the  transmit  probe,  k  as  the  propagation  vector  of 
the  scattered  plane  vave  (to  the  receive  probe),  and  r^  as  the 
coordinates  of  the  receive  probe. 

The  measurement  equation  for  near-field  bistatic  RCS  is 
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bo<'Vr2>=  V1)'  lll<k-1>W1>e 


-cIt  +tk-r, 

ljf*  4 

dL  e 


(2)  {3}  {  «) 


where 


K  =  k  e  +  k  e  (Oa) 

x  x  y  y  '  ' 

L  =  1  e  +  1  e  .  (8b) 

x  x  y  y  '  ' 

The  scalar  b^  is  a  measurable  quantity  that  is  the  conplex 
gain  factor  describing  the  transmission  path  from  the  transmit 
probe  to  the  receive  probe  via  the  scattering  object.  The 
incident  wave  is  represented  by  the  factor  {3}.  The  exponential 
{4}  represents  propagation  of  the  incident  wave  from  the 
transmitting  antenna  to  the  object,  and  the  multiplication  with 
{2}  represents  scattering  by  the  object.  The  resulting  scattered 
spectrum  of  waves  propagates  {5}  from  the  object  to  the  receive 
probe;  {1}  represents  the  interaction  of  the  propagated,  scattered 
wave  with  the  receiving  probe  itself.  Note  that  nultiple 
interactions  are  not  included  in  the  neasureTTent  equation. 

If  we  consider  the  scattering  object  as  a  transmitting 
antenna,  the  conplicated  wavefront  generated  by  it  can  be 
represented  by  a  superposition  of  plane  waves.  We  denote  that 
superposition  by  F*  (r^ ,k ) ,  vfiere  r^  is  the  location  of  the  actual 
transmit  source,  and  k  is  the  propagation  vector  of  the  scattered 


plane  waves. 
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Since  the  scattered  spectrum  is 


F1(r1»k)  = 


-cl -r , 


i1i<k’1>-Ii®<1>  e  -L' 


the  measurement  equation  above  becomes 


V) 


bo<rl'r2! 


1  = 


+ck-r 

Fi(rl>k)  e  dK‘ 


This  is  the  measurenent  equation  for  near-field  antenna 
measurement.  The  usual  manipulations  (i.e.,  tvo  sets  of 
neasur  sments  with  different  Fourier  transform,  and 

simultaneous  equation  solution)  are  enpioyed  to  determine 


(10) 


Furthermore ,  the  above  two  equat ions  are  simi lar  in  form, 
except  that  b  (r^,^)  is  a  scalar  vhile  F^r^fc)  is  a  vector.  The 
solution  techniques  of  near-field  antenna  measurement  can  be 
applied  to  solving  for  F^(r^,k),  except  that  tvn  independent 
solutions,  representing  the  e  and  e  components  of  F.(r  ,k),  must 

(j  y  1  1. 

be  found. 

In  order  to  calculate  the  scattering  matrix  it  is  necessary 
to  determine  a  second,  independent  scattered  spectrum  F^fr^.k) 
that  represents  the  scattered  spectrum  with  the  object  illuminated 
by  a  different  source  spectrum.  For  example,  with 
linearly-polarized  transmit  probe,  rotating  the  transmit  probe  by 
n/2  creates  a  different  source  spectrum.  So,  with  the  transmit 
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probe  antenna  rotated  n/ 2,  we  repeat  the  entire  series  of 


measurements  and  calculations.  For  each  transmit  probe  location 
r j  we  treasure  an  array  of  scalar  values  as  the  receive  probe  is 
moved  through  its  scan  pattern,  ft  secoiJL  array  of  scalar  values 
is  measured  with  the  receive  probe  in  another  orientation.  These 
two  arrays  of  treasured  data,  called  bQ3(r ^r^)  an£*  t3Q4^i,r2^’  are 
used  to  calculate  F^r^k). 

Referring  to  the  general  bistatic  measurement  equation  (7), 
we  have  calculated  F^r^k)  and  F^r^.k),  x-iiere 


-tl -r. 


^(rj.k)  =  Iu(k.l)-I10(l)  * 


(Ala) 


-tl -r , 


F2(rrk)  =  li1(k* x) ' Ii®( A>  e 


are  Fourier  integrals  that  tray  be  inverted  to  gixre: 


(lib) 


iii<k'1>'Ii0<1>  =  “V  I  FiCri-A)  « 

4n 


♦tl • r . 


In  (k» 1 )  ’  xi®( 1 )  =  -V  F2(rlfi) 

4n 


hIt, 


(12a) 


(12b) 


This  time  the  integrals  are  over  all  positions  of  the 
transmit  probe,  holding  the  receive  probe  position  fixed.  Since 
the  transmit  probe  moires  in  the  plane  z  =  d,  let 

_  A.  A  A 

R.  =  r,  e  +  r,  e  +  d-e  ,  (l: 

1  lx  x  ly  y  z’  x 

and  the  integrals  becorre 
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explicit  (kfl)  dependence, 


vfiich  has  solutions: 

01'  -  0‘ I 

v§10(i  10^ 

ln88  =  A 

01*  -  O* I 

w8  10*  w8 10* 

~  a 

o  i •  -  O' i 

yfllW  y8  10* 

1 11^(9  =  A 


(17) 

(19a) 

(18b) 

(18c) 
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I, 


01'  -  O'l 

ys  1®*  Vw 


11^  4 

provided  the  determinant 


(I8d; 


(19) 


is  nonzero  for  particular  (k,l).  This  oonpletes  the  calculation 
of  the  scattering  n»tnx. 

4.0  ANGULAR  FILTERING 

Since  the  probe  antennas  may  have  very  broad  patterns, 
probe-to-prube  coupling  is  an  area  of  particular  concern  in  this 
technique  of  BRCS  measurement .  However,  the  transmit  probe  is  in 
the  plane  of  the  receive  probe,  so  the  plane  wave  spectrum 
illuminating  the  receive  probe  contains  priirariiy  waves  with 

=  0.  These  plane  waves  appear  on  the  periphery  of  the  Fourier 
transform  in  k^-k^  space  and  can  be  suppressed,  provided  that  the 
probe -to -probe  coupling  has  not  swanped  the  desired  scattered 
signal  due  to  limited  receiver  dynamic  range. 

In  Figure  2  we  show  a  typical  measured  angular  spectrum.  The 
scattering  object  was  a  6- inch  diameter  sphere  illuminated  at  10 
Ghz  by  a  stationary  transmit  probe.  The  k^-k^  plane  shown  covers 
the  entire  back scattering  hemisphere,  but  the  angles  of  valid 
reconstruction  enconpass  only  the  large  double  peaks,  viiich  are 
one  coirponent  of  the  plane  wave  spectrum  scattered  from  the 
sphere.  The  broad  lower  peak  behind  the  double  peaks  comes  from 
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probe  coupling,  vfule  the  smaller  peaks  around  the  periphery  are 
due  to  room  reflections.  Note  that  the  probe  coupling  peak  has 
wrapped  around  into  the  posit ive-k^  area  of  the  calculated 
spectrum. 

5.0  DISCRETIZATION  ISSUES 

As  in  the  case  of  near-field  antenna  measurement,  it  can  be 
shown'’  that  I(k)  in  the  x-y  plane  is  effectively  hand limited,  so 
that  a  sampled-data  approach  to  measuring  I(k)  in  the  plane  z  =  d 
(a  constant)  is  sufficient  for  recons true ting  I(k).  If  the  plane 
of  measurement  is  not  in  the  reactive  near-field  of  any  source 
antennas  (i.e.,  several  wavelengths  or  more  from  the  antennas), 
sampling  I(k)  at  intervals  slightly  smaller  than  X/2  ensures  that 
evanescent  waves  are  attenuated  sufficiently  so  that  the  spatial 
Nyquist  requirement  is  met.  Finally,  spatial  filtering  and 
aliasing  theory  (and  experiment)  indicate  that  sampling  at 
intervals  larger  than  X/2  may  be  acceptable  in  certain  situations. 

Evaluation  of  integrals  of  the  form 

i*5  i 

£.k  x 

H(kx)  =  h(x)  e  M  dx  (20) 

J  —CO 

is  feasible  for  handlimited  h(x)  using  sampled  data  and  discrete 
Fourier  transform  techniques.  The  continuous  function  h(x)  is 
sampled  at  intervals  of  6x,  viiere  6x  is  chosen  sufficiently  small 
to  ensure  that  the  spectrum  of  h(x),  convolved  with  the  sampling 
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window,  is  unaliased  for  the  k  of  interest.  The  sanpled  h(x)  is 
a  f inite  sequence 

h'(n)  =  (h(x)  |  x  =  xq  +  (n-l)-5x,  n=l,...N)  (21) 

=  h* ( 1 ) ,h' (2) , . . .h' (N) 

which  is  implicitly  periodic  in  n  with  period  N.  Tte  DFT  of  h'(n) 


li 

H'(m)  =  ^  h*  (i.)t 


6  (n-i )  m-j 


,  m  =  1 ,  . . .N 


where  (for  even  N) 


ci  2n 

6k  = 
x  N6x 


m  =  1,  2, 


|  [m-(M+l)]6kx  n»  =  -4p*2,-|-^3,  ...  h) 

H(kR)  =  6x  -H‘ (tn) .  (2? 

Die  to  the  discrete  implementation,  laboratory  measurements 
and  data  analysis  yield  a  set  of  scattering  matrices 


IniW  * 


I.,  (K . ,L  . )  I,,  (K.,L  ) 

llxxv  i  j'  llxy'  i  j' 

I.,  (K.,L.)  Itl  (K.,L.) 

1 lyx v  i  j '  1 lyy '  i ’  j ' 


where  K.  and  L.  can  take  on  only  a  discrete  set  of  values  that  are 
^  J 

determined  fay  the  measurement  grid  size,  the  spatial  sampling 
interval  6x,  and  the  "rules"  of  the  DFT.  Fourier  interpolation 
can  be  used  to  increase  the  resolution  of  the  grid  of  values  of  K 


The  scan  table  has  two  movable  carriages,  each  with  one 
movable  platform.  The  carriages  and  platforms  ride  on  stainless 
steel  rails  with  recirculat ing  bai 1-bearings ,  and  are  driven  by 
Slo-Syn  stepper  motors  via  steel-and-plast ic  chains.  Step  size  is 
©.@1175  inch,  and  positioning  repeatability  is  typically  0,005 
inch.  The  table  is  equipped  with  safety  limit  switches  and 


infrared  indexing  sensors.  The  object  under  test  (OUT)  is 
supported  on  a  test  fixture  above  the  scanning  table,  and  the 
transmit  and  receive  probe  antennas  are  mounted  on  the  platforms. 
Absorber  material  is  suspended  above  the  OUT. 

The  probe  antennas  are  identical  equal -length  open-ended 
sections  of  X-band  vaveguide.  Each  probe  is  held  in  a  nounting 
bracket  that  is  bolted  to  a  single  mounting  hole  on  its  platform. 

The  probe  antennas  trove  through  predetermined  scan  patterns 
viien  measurements  are  being  made.  The  scan  patterns  are  (at  the 
moment )  non -over lapping  squares.  During  "on-the-f ly"  scanning  the 
receive  probe  mows  along  the  x-axis  while  holding  v  constant. 

The  x-  and  y-sanpling  intervals  for  both  the  receive  and  transmit 
probes  are  identical  (after  spatial  filtering  and  r^conpl ing ) . 

The  equipment  conf iguration  for  performing  near  field 
measurements  is  shown  in  Figure  4.  The  synthesized  RF  signal  is 
taken  from  the  output  port  of  the  Hewlett -Packard  8408B  Microveve 
Network  Analyzer  to  a  20-watt  TWT  amplifier,  and  then  to  the 
transmit  probe  antenna  via  senu-rigid  coaxial  cable.  Semi-rigid 
coax  also  carries  the  received  signal  from  the  receive  probe  to  a 
preamplifier  and  then  to  the  input  port  of  the  8408B. 

Measurements  are  obtained  "on  the  fly"  as  the  receive  prove 
scans  at  30  cm/ sec.  The  gain  and  phase  values  are  captured  by 
dual  sanple/hold  amplifiers  with  a  6  f/sec  window  for  digitization 


334 


by  a  12-bit/D  converter.  The  network  analyzer  bandwidth  is 

f  =  10  Khz;  assuming  an  effect ive  integration  t me  of  about 

3/f.  ,  the  data  "smear"  length  is  about  0.  L  inn. 

bw 

Measured  gain  and  phase  are  con^rted  to  I  and  Q  conponents 
and  stored  to  hard  disk.  Data  acquisition  is  controlled  by  a 
conputer  program  that  runs  in  the  IBM-PC  computer. 

7.0  OVER-SAMPLING 

Under  ideal  conditions  the  gain  and  phase  need  be  sampled 
only  at  A/2  internals  (or  slightly  more  often  if  one  is  near  the 
reactive  near-field  of  the  test  object).  Ue  oversample  by  a 
factor  of  ten  because:  1)  noise  generated  in  the  electronics  for 
gam  and  phase  detection  smears  out  the  spectral  content  of  the 
signal  being  measured;  2)  a  general  rule  of  thumb  in  digitizing 
and  processing  noisy  signals  is  that  one  should  digitize  at  five 
to  ten  times  the  Nyquist  rate;  3)  filtering  techniques  can  be 
used  to  reduce  the  (uncorrelated )  noise  on  the  signal; 

4)  nonlinear  filtering  techniques  can  be  used  to  detect  and 
correct  invalid  phase  measurements  caused  by  phase  wrap-around, 
and  5)  room  reflections  produce  spectral  components  with  spatial 
frequencies  very  near  the  spatial  Nyquist  frequency.  The  measured 
noise  will  be  unrorrelated  from  sample  to  sample  if  digitizing  is 
performed  at  a  rate  that  is  at  least  an  order  of  magnitude  slower 
than  tiie  reciprocal  of  the  IF  bandwidth. 
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The  data  are  sarrpled  at  \/20  intervals  dur  1  ng  scanning.  The 
algorithm  to  detect  glitches  that  occur  in  the  treasured  phase 
values  (due  to  phase  virap-around  in  the  gain/phase  receiver)  uses 

differentiation  to  enhance  noise  and  attenuate  low- frequency 

6  7 

deterministic  signal  content  ,  and  Tchebychef f ’ s  Inequality  to 

detect  unusually  large  iVj i %e  peaks  uhat  are  due  to  the  phase 

glitches.  Tchebychef f ’ s  Inequality  relates  the  probability  of 

deviation  of  a  sample  of  a  stochastic  process  from  its  expected 

value,  to  the  variance  of  the  process,  as 

P  {  )  X— >u )  >  no}  <  n  2  (3©) 

.jhere 

f*  5  E{X)  (31) 

o  =  e|(X-/i)2J  (32) 

are  the  expected  mean  and  variances,  respectively.  For  example, 
if  X  is  the  I  signal’s  derivative,  and  n=3,  then  X  will  deviate 
from  its  mean  by  an  arount  3o  with  probability  no  larger  than 
0.111,  irrespect  i<re  of  the  probability  density  function  of  the 
stochastic  process.  With  appropriate  choice  of  n,  this  technique 
detects  phase  glitches  reliably  without  triggering  on  noise 
components  that  are  indeed  part  of  the  receiver’s  noisy  output 
signal . 

The  detected  phase  glitches  in  I  and  Q  are  patched  over  by 
interpolat ion  and  a  low- pass  spatial  filter  is  applied  to  I  and  Q. 
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While  averaging  groups  of  ten  adjacent  samples  is  a  straight 
forward  approach  to  filtering,  it  distorts  the  higher  spatial 
frequency  components  of  the  desired  signal.  We  use  a  21-term 
linear-phase  low-pass  f inite-inpulse-response  filter  which  is 
applied  to  the  over-sampled  rawr  data  to  improve  the  signal  to 
noise  ratio.  The  filter  response  is  down  ©.?  db  at  spatial 
frequencies  of  i/A,  and  5®  db  down  at  4/X.  The  filtered  data  are 
resanpled  at  \'2  and  stored  to  disk.  Spatial  filtering  inposes  no 
time  performance  penalty  since  the  computer  is  idle  during  the 
scan  re - 1  race  i nt  er va 1 . 

By  means  of  power  spectral  analysis  of  the  raw’  and  filtered 
data  sets  we  determined  that  the  I  and  Q  signals  are  dovn  at  least 
25  db  at  !/>.,  and  that  phase  glitch  correction  and  spatial 
filtering  lower  the  effect ive  noise  floor  by  typically  1®  to  15 
db. 

8.0  VERIFICATION  APPROACH 

Since  the  experimental ly-determined  scattering  matrices  are 
defined  at  discrete  values  of  K  and  L,  it  is  convenient  to 
construct  one’s  theoretical  models  of  scattering  such  that  the 
precise  values  of  K  and  L  calculated  by  the  analysis  software  can 
be  automatically  plugged  in  to  yield  the  theoretical  scattering 
values.  From  an  algorithmic  standpoint,  this  corresponds  to 
implementing  the  theoretical  or  numerical  model  of  scattering  as  a 


subrout ine  that  has  as  input  the  values  of  K  and  L  for  which  a 
theoretical  prediction  is  needed. 

Vfe  are  taking  a  three-part  approach  to  verifying  the 
equipment  and  theory.  Ir  the  first  part,  we  per f or rred  a  standard 
near-field  antenna  Treasure vnent  on  a  knovn  antenna,  processed  the 
data  using  analytic  probe  pattern  for  probe  correction,  and 
compared  the  resulting  pattern  data  to  known  data  for  the  antenna. 
The  result  of  this  step  was  verification  of  the  near-field  antenna 
measurement  capability  {equipment,  data  acquisition  and  analysis 
software).  Since  the  first  part  of  the  BRCS  data  analysis  is 
essentially  a  near-field  antenna  solution,  this  provided  a  partial 
validation  of  the  SICS  software. 

In  the  second  part  we  performed  near-field  antenna 

measurement  on  a  known  target  (conducting  sphere)  with  known 

qua s l -p lane- wave  illumination  in  the  E  and  E,  planes.  The  usual 

a  9 

near-field  antenna  data  analysis  was  performed  using  a 
high-density  grid,  and  *he  results  compared  with  predictions.  The 
result  of  this  part  was  verification  of  the  prediction  and 
comparison  sof t ware . 

The  third  part  consists  of  performing  the  full-up  bistatic 
scan  with  both  transmit  and  receive  probe  antennas  in  motion.  The 
measured  data  will  be  processed  as  described  above  to  yield  the 
scattering  matrix  for  the  object.  The  results  of  this  part  will 
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be  verification  of  the  bistatic  analysis  software  and  an 
evaluation  of  the  feasibility  of  the  PW3M  bistatic  RCS  algorithm. 
Both  separate  and  overlapping  scan  patterns  will  be  evaluated. 
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Figure  1  —  Coordinate  System 
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Figure  2  —  Typical  Measured  Angular  Spectrum 
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Figure  4  —  Apparatus 
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ABSTRACT 

This  paper  discusses  a  criterion  useful  in  the  design  of  a 
two-wire  simulator  in  order  to  optimize  the  uniformity  of  the 
incident  FM  fields  within  the  volume  occupied  by  a  test  object. 
It  is  shown  that  the  separation  between  the  two  wires  2a  and  their 
height  above  the  ground  b  can  be  chosen  in  such  a  way  to  minimize 
a  functional  describing  the  deviation  from  uniformity,  in  the 
sense  of  the  root-mean-square  value,  of  the  principal  component  of 
the  electric  field  in  a  cross  section  of  the  volume  occupied  by 
the  test  object.  Measurements  performed  on  an  actual  two-wire 
illuminator  verify  that,  when  this  criterion  is  implemented  in  the 
design,  the  field  uniformity  is  satisfactory. 
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1.  INTRODUCTION 


The  purpose  of  the  analysis  described  in  this  paper  is  to 
support  the  design  of  an  illuminator  for  EMP  ( El  ectromagnetic 
Pulse)  hardness  surveillance.  The  antenna,  shown  in  Figure  1, 
consists  of  two  wires  in  a  configuration  similar  to  a  rhombus 
connected  i.u  a  generator  aL  one  vertex  and  to  a  matched  termina¬ 
tion  load  at  the  other  vertex.  At  the  section  PP'  the  wires  are 
supported  by  a  dielectric  cord  stretched  between  two  poles.  The 
mechanical  arrangement  allows  for  the  wire  separation  and  height 
above  the  ground  to  be  varied  within  an  ample  range,  thus  offering 
design  flexibility.  The  portion  of  the  system  on  the  left  of 
section  PP1  will  be  referred  to  as  the  launch  region,  whereas  that 
on  the  right  is  the  termination  region.  The  otject  to  be  tested 
occupies  the  volume  indicated  by  V,  which  will  be  referred  to  as 
the  working  volume.  The  generator  can  drive  the  wires  at  either 
the  same  voltage  (common  mode)  or  at  voltages  of  opposite  polarity 
(differential  mode). 

When  the  structure  rests  on  a  perfectly  conducting  ground  and 
has  a  matched  termination,  it  has  been  shown  [1],  [2],  for  the 
common  mode  excitation  that,  at  low  frequencies,  the  currents 
along  the  wires  are  well  approximated  by  travelling  waves  and  that 
the  system  behaves  like  a  two-wire  transmission  line  with  a  spher¬ 
ical  wavefront  operating  in  the  TEM  mode.  Furthermore,  since  in 
the  present  case  the  angle  n  indicated  in  Figure  1  is  maintained 
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small  in  the  launch  region,  the  portion  between  wires  and  ground 
of  the  surface  of  the  sphere  (wavefront)  which  has  the  source  as 
origin  and  z  as  radius  (see  Fig.  1)  can,  with  a  good  approxima¬ 
tion,  be  taken  to  be  its  projection  on  the  plane  x-y. 

As  the  frequency  is  increased  the  travelling-wave  transmis¬ 
sion  line  model  is  not  truly  valid  anymore  and  both  the  current 
and  the  fields  exhibit  a  standing  wave  behavior.  Moreover,  longi¬ 
tudinal  components  of  the  electromagnetic  fields  appear  in  addi¬ 
tion  to  transverse  components.  However,  for  most  of  the  frequency 
range  of  interest  to  us,  on  each  transverse  plane,  the  formal 
distribution  of  the  dominant  component  of  the  fields  can  still  be 
reasonably  approximated  by  those  of  a  transmission  line.  The  EM 
fields  of  a  two-paral  1  el -wi  re  transmission  line  above  a  perfectly 
conducting  plane  operating  in  the  TEM  mode  in  any  plane  transverse 
to  the  z-axis  are  well  known  from  static  analysis.  Such  distribu¬ 
tion  is  sketched  in  Figure  ?  for  both  the  common  mode  and  the 
differential  mode.  It  is  observed  that,  in  a  region  D  between  the 
wires  and  the  yround  which  is  a  cross  section  of  the  working 
volume,  for  the  common  mode  excitation  the  electric  field  is 
predominantly  vertical,  whereas  for  the  differential  mode  it  is 
mostly  horizontal.  When  performing  EMP  hardness  surveillance 
tests  one  is  concerned  with  illuminating  a  test  object  with  elec¬ 
tromagnetic'  fields  which  are  as  close  as  possible  to  those  of  a 
plane  wave,  i  .e. ,  uniform  in  any  transverse  plane.  Therefore,  the 


main  goal  of  this  analysis  was  to  determine  the  geometrical  con¬ 
figuration  of  the  wires,  i .e. ,  separation  2a  and  height  b  above 
the  ground  (see  Fig.  2)  which  give  the  most  uniform  field  in  a 
prescribed  region  D,  for  both  common  mode  and  differential  mode 
operations. 

2.  MATHEMATICAL  FORMULATION 

The  transverse  EM  fields  of  a  transmission  line  operating  in 
the  TEM  mode  can  be  determined  by  solving  the  static  problem  of 
two  infinitely  long,  uniformly  charged,  parallel  wires  above  a 
perfectly  conducting  ground  plane.  This  two-dimensional  problem 
is  described  hy  the  Cartesian  coordinate  system  presented  in 
Figure  2.  The  distribution  of  the  fields  everywhere  in  the  x-y 
plane  is  easily  obtained  by  applying  the  method  of  images  and 
superposition.  When  the  charge  has  the  same  sign  on  both  wires 
(corresponding  to  the  common  mode)  the  electric  field  in  the 
region  D  is  predominantly  directed  along  the  y-axis.  Such  domin¬ 
ant  component  is  given  by 


E  =  -J—  ! 

y  2irr 


y  -  b 
1 


7 


y  -  b 

T 


o  (x-a)  +  (y-b)  (x+a)  +  (y-b 


_ y  +  b _  _ y  +  b _ 

(x-a)'  +  (y+b)  (x+a)^  +  (y+b) 


(1) 


where  q  is  the  wire  charge  per  unit  length  and  zQ  is  the  vacuum 


electric  permittivity.  On  the  contrary,  when  the  wires  have 
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opposite  charge  (corresponding  to  the  differential  mode)  the 
electric  field  in  D  is  predominantly  directed  along  the  x-axis. 
Such  dominant  component  is 


F  _  q _  f  x  +  a  x  -  a 

x  K  (x+a)2  +  (y-b)2  (x-a)2  +  ( y-b) 2 

.  -  .  x  _  +  .  x  -  a  i  f2x 

(x+a)2  +  (y+b)2  (x-a)2  +  (y+b)2 

Within  a  prescribed  region  D  the  uniformity  of  the  field  depends 
on  the  choice  of  a  and  b.  As  a  criterion  to  quantify  the  extent 
of  uniformity  we  introduce  a  functional  F  (E^,  a,  b,  D)  [3],  which 
gives  the  relative  deviation,  i  .e. ,  the  2-norm  error,  of  any  of 
the  relevant  components  of  the  field  with  respect  to  the  average 

value  taken  in  D.  In  its  general  form  F  (E.  ,  a,  b,  D)  is  defined 

<5 

as 


F(Ef,a,b,D) 


=  li 


r  (a  ,b  ,£ )  -  E^ve(a,b)l2dd1/2/E^ave(a,b) 


(3) 


where  F.  is  either  Eu  for  the  common  mode  or  E„  for  the  differen- 

rS  J  * 

tial  mode  and  r  is  a  variable  describing  the  domain  D. 
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3  yp 

Moreover,  E  is  defined  as 

Eaye  1  f  p  (a,b,0  d,  (4) 

6  U  D  6 

Since  we  are  concerned  here  with  test  objects  such  as  aircraft 

whose  cross  sections  exhibit  a  maximum  linear  dimension  which  is 
greater  than  the  others,  D  can  be  taken  to  be  a  straight  segment 
either  lying  parallel  to  the  x-axis  (representing  the  wings)  or 
along  the  y-axis  (representing  the  cross  section  of  the  fuselage 
and  tail)  at  a  certain  height  h  above  the  ground.  This  choice 
allows  us  to  carry  out  the  integrals  in  Eq.  (3)  and  Eq.  (4)  analy¬ 
tically.  Once  F  is  constructed,  the  maximum  uniformity  is  achiev¬ 
ed  by  choosing  a  and  b  such  that  F  is  minimum  in  D.  Moreover, 

since  a  and  b  are  independent  parameters  it  might  be  useful  to 
introduce  an  additional  constraint  based,  for  instance,  on  the 
behavior  of  the  field  or  its  derivatives  at  one  point,  to  esta¬ 
blish  a  relationship  between  a  and  b,  thus  eliminating  one  para¬ 
meter.  For  the  case  of  the  common  mode  excitation  one  such  rela¬ 
tionship  exists  T 43  that  provides  a  good  uniformity  in  a  local 
region  around  the  center  of  the  structure,  coincident  with  the 
origin  of  the  coordinate  system  (from  Fig.  2),  where  the  field  is 
strong.  Such  relationship  is  a/b  =  1//J  and  it  has  been  incorpor¬ 
ated  in  our  analysis.  For  the  differential  mode  excitation  a 
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correspondent  relationship  was  not  found  because  in  this  case  the 
origin  is  a  point  when  the  field  goes  to  zero. 

Once  the  integrations  in  Eq.  (3)  and  Eq.  (4)  are  carried  out 
one  is  left  with  transcendental  functions  of  a  and  b  which  are 
then  minimized  graphically. 

For  the  sake  of  brevity  the  explicit  expressions  are  not  reported 
here. 

Figure  3  presents  the  resulting  2-norm  error  over  the  wings 
for  the  common  mode  excitation,  as  a  function  of  a/h.  The  rela¬ 
tionship  a/h  =  1//3  is  here  used.  It  is  seen  that  the  error 
decreases  rapidly  and  for  a/h  >  5  it  is  bounded  by  5  percent. 
From  this  curve,  for  any  prescribed  maximum  deviation,  one  can 
immediately  determine  a  range  of  possible  values  for  a,  and  conse¬ 
quently  for  h,  which  satisfies  the  requirement. 

Figure  4  illustrates  the  2-norm  error,  again  over  the  wings, 
for  the  differential  mode  excitation.  In  this  case,  one  is  con¬ 
cerned  with  a  family  of  curves  with  the  variable  b/h  on  the  ab¬ 
scissa  axis  while  a/h  acts  as  a  parameter,  assuming  discrete 
values.  A  dual  representation,  with  the  roles  of  a  and  b  inter¬ 
changed,  provides  exactly  the  same  information.  Again  it  is  noted 
that  the  error  decreases  when  a/h  and  b/h  increase.  Also,  for  a/h 
and  b/h  large,  which  represent  a  realistic  situation,  the  minima 
of  the  curves  are  fairly  broad,  emphasizing  that  the  choice  of  a 
and  b  is  not  critical  within  a  certain  range  of  values.  Charts 
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Tike  those  of  Figure  3  and  4  are  directly  useful  to  the  designer 
in  choosing  the  optimum  a  and  b  within  a  range  of  possible  values, 
compatibly  with  other  constraints,  such  as  test  site  area,  pole 
length,  etc. 

3.  IMPEDANCE  CALCULATIONS 

For  good  performance,  illuminators  of  this  type  are  usually 
terminated  into  a  resistive  load  of  value  equal  to  the  character¬ 
istic  impedance  Zc  =  V/I  of  the  transmission  line  at  the  ending 
section.  Using  again  the  static  model  for  two  parallel  wires  to 
calculate  the  voltage  and  after  dividing  by  the  current,  it  is 
obtained 


V  = 


4^ 


4n  f 


1  +  /l+(R/b)2  -  R/b-,  x  1 


1  -  /I+(R/b)2  +  R/b 


■1  [1  +  (b/a)  ] }  (5) 


for  the  common  mode  and 


Z"  =  ^  h. n  r1  *  /T^R/a^---R-/Al  1  ln  fl  +  ( a/b)2] }  (6) 

c  v  1  -  /l+(R/a) i  +  R/a 

for  the  differential  mode,  where  R  is  the  radius  of  the  wire 
and  n  =  377  a  is  the  free-space  intrinsic  impedance.  It  is 
pointed  out  that  in  both  cases  the  characteristic  impedance  con¬ 
sists  of  the  sum  of  two  terms;  one,  depending  either  on  the  ratio 
R/b  (common  mode)  or  R/a  (differential  mode)  is  the  dominant 
contribution,  while  the  other,  depending  on  the  ratio  a/b,  pro- 
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impedance  at  PP * .  This  can  be  of  great  help  in  reducing  reflec¬ 
tions.  Unfortunately,  this  correction  cannot  be  applied  to  the 
launch  region  because  here  d/dz  {Zc>  is  positive  and  one  would 
have  to  "sdbstract"  a  resistance  d/dz  {Zc>  dz  for  every  portion  dz 
of  wires. 

It  is  pointed  out  that  the  reflections  and  scattering  origin¬ 
ating  at  the  wire  bends  are  not  affected  by  this  measure  and  hence 
they  are  still  present. 

One  might  try  to  mitigate  this  effect  by  making  a  smooth  transi¬ 
tion  around  the  poles  instead  of  having  a  sharp  bend. 

4.  CONCLUSIONS 

A  criterion  is  presented  to  help  in  selecting  the  wire  separ¬ 
ation  and  height  above  the  ground  in  the  design  of  a  rhombic 
illuminator,  such  that  the  most  uniform  fields  in  a  prescribed 
working  volume  are  obtained.  In  addition  some  impedance  consider¬ 
ations  which  help  to  reduce  the  reflections  in  the  working  volume 
are  also  given.  The  above  criterion  and  impedance  considerations 
were  applied  in  the  design  of  a  low  level  CW  illuminator.  Unlike 
the  configuration  in  the  model,  the  real  system  was  built  without 
a  conductive  ground  plane  in  order  to  get  higher  fields  near  the 
ground  for  the  differential  mode  excitation.  Results  of  a  prelim¬ 
inary  field  mapping  for  the  common  mode  excitation  show  that  the 
uniformity  is  satisfactory  and  that  the  distribution  of  the  prin¬ 
cipal  component  of  the  electric  field  in  any  cross  section  of  the 


working  volume  closely  follows  that  of  the  model  used  in  the 
analysis  (Eq.  1),  with  differences  of  less  than  a  factor  of  two. 
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Abstract 

The  ongoing  development  of  a  computer-aided  design  tool  for  the 
comprehensive  analysis  of  antenna  systems  and  subsystems  is  described.  Utilizing 
custom  Automated  Network  Analyzer  software,  arbitrary  N-port  devices  are  fully 
characterized  in  terms  of  scattering  coefficients.  This  data  is  in  turn  modeled  as  a 
function  of  frequency  and  physical  variables  in  order  to  provide  a  continuous 
parametric  representation  with  up  to  a  50:1  reduction  in  required  storage  space. 
This  parametric  representation  allows  for  a  statistical  investigation  of  antenna 
system  and  subsystem  behavior  in  terms  of  user-specified  mechanical  and  electrical 
tolerances.  Based  on  an  accrued  library  of  scattering  matrix  characterizations  of 
individual  antenna  components,  the  resultant  network  parameters  (VSWR,  Isolation, 
Insertion  Loss,  and  Phase),  far-field  patterns,  and  RCS  characteristics  of  ensemble 
antenna  systems  and  subsystems  can  be  predicted  based  on  simple  linear  algebra 
techniques.  The  analytical  theory  and  mechanics  of  this  approach  are  discussed  along 
with  numerical  results  for  a  typical  antenna  subsystem  application. 
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1. 


Introduction 


In  an  environment  of  increasingly  stringent  performance  and  observability 
requirements  across  ever-broadening  frequency  bandwidths,  today's  antenna  designer  must 
reconsider  the  soundness  of  heretofore  common  analysis  assumptions  and  approximations. 
Specifically,  higher-order  reflections  and  finite  isolation  effects  between  adjacent  and  remote 
components  can  no  longer  be  dismissed  as  negligible,  the  finite  accuracy  of  theoretical 
models  becomes  a  point  of  increasing  concern,  the  imperfections  of  "perfect"  loads  take  on 
new  significance.  Likewise,  the  prudent  specification  of  internal  mechanical  and  electrical 
tolerances  becomes  increasingly  critical  in  optimizing  the  balance  between  performance  and 
producibility.  In  response  to  these  new  challenges,  we  are  currently  developing  a  unified 
approach  to  the  analysis  of  antenna  subsystem  and  system  behavior  based  on  the 
measurement  and  utilization  of  parameterized  scattering  coefficient  data. 

Scattering  coefficients  were  chosen  over  other  characterizations  (Impedance, 
Admittance,  etc.),  because  of  their  straightforward  means  of  measurement  and  their 
correspondence  to  real  physical  attributes.  Figure  1  illustrates  an  overall  flow  diagram  of  our 
approach  which  can  be  conveniently  subdivided  into  two  major  subtasks. 

First,  as  a  resource  for  the  second  task,  we  create  and  augment  the  scattering 
parameter  library.  This  involves  the  scattering  coefficient  characterization  of  a  variety  of 
common  devices  and  junctions  in  terms  of  frequency  and  physical  variables  (iris  depths, 
septum  offsets,  stripline  widths,  etc.)  over  suitable  ranges.  This  raw  data  is  first  collected 
through  the  utilization  of  an  automated  test  procedure,  transformed  and  processed  in  the 
spatial  domain,  parameterized  in  terms  of  physical  variables,  and  finally  stored  in  condensed 
form  as  a  member  of  the  scattering  parameter  library. 

The  second  major  subtask  is  the  computation  of  the  composite  scattering  coefficients 
(and  related  characteristics)  for  user-specified  networks  comprised  of  components  whose 
individual  scattering  coefficients  have  been  previously  stored  (in  condensed  form)  in  the 
scattering  parameter  library.  These  computations  are  performed  over  a  user-specified  range 
of  frequencies  for  a  user-specified  set  of  physical  variable  settings  and  tolerances.  An 
executive  program  manages  overall  analysis  operations  in  conjunction  with  a  circuit  analysis 
"engine"  optimized  for  the  performance  of  repetitive  linear  algebra  operations. 
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DEVICE 


FIGURE  1.  OVERALL  FLOW  DIAGRAM,  THE  UNIFIED  ANALYSIS  APPROACH 
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2. 


Creation  and  Augmentation  of  the  Scattering  Library 


2.1  Scattering  Coefficient  Measurement 

Figure  2  illustrates  the  fully-developed  automated  test  procedure  in  flow  diagram 
form.  This  procedure  is  currently  software-implemented  on  an  HP8409C  Automated  Network 
Analyzer  as  depicted  in  Figure  3,  Data  is  temporarily  stored  on  floppy  disk  at  the  HP9000 
terminal  prior  to  transmission  to  a  mainframe  (TSO)  computer  via  a  modem  switch. 

The  automated  test  procedure  commences  with  a  series  of  prompts,  at  which  point 
the  user  provides  a  description  of  the  device  under  test  (device  function,  reciprocity,  number 
of  ports,  dataset  name.)  In  addition,  interface  specifications,  interface-to-junction  distances, 
and  names  are  specified  by  the  user  for  each  device  port.  The  prompting  sequence 
continues  with  user-specification  of  lower  limit,  upper  limit,  and  cut-off  frequencies  plus  the 
desired  number  of  frequency  samples  (restricted  to  an  integer  power  of  2  due  to  the 
implementation  of  Fast  Fourier  Transform  (FFT)  in  later  processing.)  Based  on  this 
information,  frequency  samples  are  selected  such  that  they  are  evenly-spaced  in  phase 
constant  p  in  order  to  provide  optimal  resolution  when  transformed  to  the  spatial  domain  (see 
Section  2.2.)  Finally,  the  user  is  prompted  for  the  names,  limits,  and  increments  of  all  physical 
variables  to  be  associated  with  the  device. 

Based  on  the  information  provided  by  the  user,  an  optimized  test  sequence  is 
formulated  and  presented  to  the  user  for  approval.  This  tabular  test  plan  specifies  the 
minimum  number  of  required  permutations  of  port  connections  and  physical  variable 
combinations,  minimized  on  the  basis  of  device  reciprocity  and  network  analyzer  capabilities. 
In  addition,  data  storage  space  is  allocated  and  initialized  with  all  user-specified  device 
information. 

Upon  user  approval,  a  calibration  procedure  is  performed  and  the  test  sequence  is 
initiated.  All  required  hardware  disconnections  and  connections  are  explicity  prompted  and 
confirmed  step-by-step. 

All  raw  measured  data  is  calibrated  at  the  completion  of  the  test  sequence,  then 
formatted  prior  to  transfer  to  the  mainframe  computer.  Transfer  is  accomplished  via  a  modem 
switch  with  all  data  echoed  for  confirmation  of  error-free  transmission. 


364 


REFL  TRANS. 


DEVICE 
UNDER  TEST 


FIGURE  3.  TEST  SET-UP,  AUTOMATED  TEST  PROCEDURE 
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Data  Reduction  in  the  Spatial  Domain 


The  overall  methodology  of  this  facet  of  the  modeling  process  is  illustrated  in  Figure 
4.  However,  before  discussing  this  approach  in  detail,  it  seems  appropriate  to  justify  our 
methods  based  on  some  theoretical  and  physical  reasoning. 

Calibrated  data  transferred  from  the  network  analyzer  to  the  mainframe  computer  is  in 
the  form  of  N-frequency  samples  (for  each  scattering  coefficient),  evenly-spaced  in  phase 
constant  (5.  Note  that  we  refer  to  p  as  the  "frequency"  domain  because  of  the  simple  one-to- 
one  relationship  between  p  and  f  (directly  proportional  for  fc  =  0.) 

The  frequency  domain,  while  the  most  direct,  is  not  the  ideal  domain  for  modeling  the 
physical  variable  dependence  of  the  individual  scattering  coefficients.  Instead,  transformation 
to  the  spatial  domain  is  desirable  in  order  to  separate  frequency  and  physical  dependencies 
and  thereby  simplify  the  modeling  procedure.  To  see  this,  assume  that  any  particular 
scattering  coefficient  can  be  expressed  "exactly"  as: 

oo 

Sm„(p.u)=X  rk<Pi)e",|,lk  0) 

krr-oo 

where  r^(P,&)  is  a  discrete  complex  reflection  or  transmission  coefficient  for  a  discrete 
frequency-(P)  and  physical  variable-(&)  dependent  discontinuity  located  a  distance  1^  from  a 
reference  plane.  These  discontinuities  may  be  real  or  fictitious  (higher-order  reflection  terms.) 
It  is  our  desire  to  solve  for  the  reflection/transmissmn  coefficients  T|<(P,2L)  at  each  of  the 
discontinuity  positions  1^  given  Smn(P,£).  However,  two  obstacles  block  our  immediate 
progress  in  that  direction. 

First,  the  frequency  dependence  of  the  reflection/transmission  coefficients  appears 
to  preclude  the  direct  application  of  a  Fourier  Transform  between  P -space  (frequency  domain) 
and  l-space  (spatial  domain),  since  Fourier  coefficients  must  be  independent  of  the  transform 
variable.  Second,  the  infinite  nature  of  the  summation  would  indicate  that  an  infinite  number 
of  frequency  samples  of  the  scattering  coefficient  would  be  required.  As  a  solution  we 
propose  now  (and  validate  later)  the  following  assumptions. 

1.  Assume  that  the  frequency  dependence  of  each  discontinuity  T^p.i)  is  weak 
compared  to  the  frequency  dependence  of  the  composite  scattering  coefficient  Smn(p,x)  and 
can  therefore  be  modeled  by  a  frequency-independent  term  F^i)  augmented  by  a  "few" 
fictitious  frequency-independent  "sideband"  terms  rr(i)  adjacent  to  rk(*). 
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FIGURE  5.  DATA  REDUCTION  IN  THE  SPATIAL  DOMAIN 
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2.  Assume  that  the  r k(20  terms  diminish  with  increasing  k  sufficiently  quickly,  such 
that  the  series  can  be  truncated  at  N  terms. 

With  these  assumptions  applied,  (1)  can  be  rewritten  as: 

N-1 

smn  (Pi.2S)  =  X  rk  e  H  [1=0.1 .2 . (N-1 )]  (2) 

k=0 


which  is  in  a  form  ideally  suited  for  application  of  the  Fast  Fourier  Transform  (FFT)  in  order  to 
solve  for  the  N  reflection/transmission  coefficients  Tk(*)  at  each  of  the  discontinuity  positions 
l|<  given  N  frequency  samples  of  the  scattering  coefficients  Smn(Pj,x). 

Given  the  N  discrete  evenly-spaced  "frequency"  samples  Pj  ,  we  can  solve  for  the  N 
discrete  evenly-spaced  discontinuity  positions  1^.  These  are  easily  shown  to  be: 


=  (N-1 )  2  7t 

k  N  (Phigh~Plow) 


[k=0,1 ,2 . (N-1)]  (3) 


where  Phigh  and  Plow  are  ,he  phase  constants  associated  with  the  highest  and  lowest  sample 
frequencies. 

Note  that  the  discrete  sampled  nature  of  Smn(Pj^)  and  the  finite  summation  limit  (N-1) 
appearing  in  (2)  imply  a  minimum  resolution  Al  and  maximum  range  lmax  in  determining 
discontinuity  locations.  These  limits  are  easily  verified  to  be 

■  =  (N-1)  2k 

N  (Phigh-Plow)  *  } 

lmax  =  (N-1  )AI  (5) 


Thus,  with  the  application  of  FFT,  we  transform  the  scattering  coefficient  data  from  the 
frequeno*4jdomain  to  the  spatial  domain,  thereby  removing  direct  frequency  dependence  and 
reducing  the  modeling  problem  to  that  of  characterizing  the  physical  variable  dependence  of 
each  reflection/transmission  (Fourier)  coefficient.  This  subject  will  be  covered  in  detail  in 
Section  2.3. 

Note  that  the  discontinuities  in  a  real  physical  device  are  usually  somewhat  localized, 
therefore  it's  not  at  all  surprising  to  discover  that  the  total  spectral  energy  (defined  as  the  sum 
of  the  squares  of  the  Fourier  coefficient  magnitudes)  is  highly  localized  as  well.  We  therefore 
can  (and  do)  exploit  this  fact  in  order  to  dramatically  reduce  the  number  of  coefficients  that 
must  be  modeled  and  stored  by  selectively  retaining  only  dominant  coefficients. 
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This  advantage  will  be  discussed  towards  the  end  of  this  section.  Additionally,  the  physical 
equivalence  between  the  Fourier  coefficients  and  the  actual  discontinuity  locations  and 
reflection  coefficients  implies  a  reasonable  expectation  that  the  physical  variable  (*) 
dependence  of  the  Fourier  coefficients  will  be  well-behaved. 

With  this  background  in  mind,  we  can  now  restate  our  spatial  domain  modeling 
approach  as  follows:  Optimally  pre-  and  post-process  the  scattering  coefficient 
data  In  such  a  way  as  to  retain  a  minimum  number  of  Fourier  coefficients 
whose  physical  variable  dependence  is  maximally  benign  while  allowing 
reconstruction  of  the  scattering  coefficients  In  the  frequency  domain  to  an 
acceptable  level  of  precision.  This  procedure  is  illustrated  pictorially  in  Figure  5.  The 
mechanics  of  the  approach  are  quite  involved,  but  include  some  basic  “building  block" 
processes. 

Center-referencing  is  performed  on  the  scattering  coefficient  data  in  the 
frequency  domain  for  the  purpose  of  removing  excessive  complex  exponential  "rotation" 
through  the  artificial  removal  of  idealized  transmission  line  lengths.  Figure  G  illustrates  this 
process  and  shows  the  resultant  translation  toward  the  reference  plane  in  the  spatial  domain. 

Unmasking  is  a  process  by  which  internal  reflection  coefficients,  attenuated  by  the 
presence  of  large  discontinuities  near  the  reference  plane,  are  restored  to  their  isolated 
magnitude  levels.  This  is  accomplished  using  a  coefficient-to-coefficient  recursion  relation 
and  has  the  effect  of  separating  the  physical  variable  dependence  of  specific  coefficients, 
thereby  improving  overall  modeling  accuracy.  Figure  7a  illustrates  the  transformed  reflection 
data  for  a  two-stub  waveguide  tuner  (pictured  in  Figure  17)  with  both  stubs  protruding  halfway 
into  the  waveguide.  Note  the  clearly  dominant  Fourier  coefficients  associated  with  the  two 
capacitive  stubs  and  the  apparent  attenuation  of  the  second  dominant  coefficient  despite  the 
equal  penetration  of  the  two  stubs.  Figure  7b  illustrates  the  spatial  domain  after  application  of 
the  unmasking  process.  Note  the  amplification  of  the  second  coefficient  and  the  (now 
apparent)  pseudo-discontinuity  corresponding  to  the  second-order  reflection  contribution. 

Focusing  is  a  method  by  which  small  idealized  transmission  line  lengths  are  added 
or  removed  in  the  frequency  domain  in  order  to  enhance  spectral  energy  concentration  in  the 
spatial  domain,  thereby  minimizing  the  number  of  Fourier  coefficients  which  must  be  modeled 
and  retained.  Optimal  focusing  is  accomplished  when  the  actual  physical  location  of  the 
dominant  discontinuity  in  the  spatial  domain  falls  directly  on  a  discrete  1^  value.  The  required 
focusing  shift  range  is  therefore  restricted  to  the  width  of  one  resolution  cell  Al.  Figure  8a 
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DISCONTINUITY  #1  (XI) 


FIGURE  7a.  BEFORE  UNMASKING, 
SPATIAL  DOMAIN 


FIGURE  7b.  AFTER  UNMASKING, 
SPATIAL  DOMAIN 


FIGURE  8a.  BEFORE  FOCUSING, 
SPATIAL  DOMAIN 


FIGURE  8b.  AFTER  FOCUSING, 
SPATIAL  DOMAIN 


illustrates  the  spatial  domain  for  the  two-stub  tuner  before  focusing.  Note  the  enhanced 
swjwutAi  sneig,  c  jr.ccm"'  t;cn  afie  j  r -  ■  > > i ■.*  <  » *  f  f^vre  Bb. 

The  reason  for  the  improvement  due  to  focusing  is  best  understood  by  considering 
the  case  of  a  single  theoretical  discontinuity  of  unit  reflection  coefficient  magnitude,  located 
precisely  at  the  I  =  I,  position.  As  illustrated  in  Figure  9,  the  frequency  response  for  this 
discontinuity,  when  transformed  to  the  spatial  domain,  yields  a  precise  unit  impulse  function  at 
the  k  =  i  position.  Now  consider  the  same  discontinuity  shifted  to  the  position  I  =  lj  +  (Al/2). 
Figure  10  illustrates  the  spatial  domain  for  this  case,  revealing  equal  and  opposite  impulses  of 
magnitude  (2/rc)  located  at  the  adjacent  positions  k  =  i  and  k  =  (i+1).  A  significant  sideband 
structure  is  now  evident,  requiring  the  retention  of  numerous  coefficients  in  order  to 
adequately  model  this  discontinuity,  and  hence  the  value  of  focusing  becomes  apparent. 

We  will  now  validate  our  previous  assumption  that  a  mildly  frequency-dependent 
discontinuity  can  be  modeled  by  a  frequency-independent  discontinuity  together  with  an 
adjacent  sideband  structure.  Figure  11  illustrates  the  spatial  domain  for  a  frequency- 
independent  discontinuity  located  at  the  physical  position  I  =  lj.  Figure  12  illustrates  the 
normalized  spatial  domain  for  a  linearly  frequency-dependent  discontinuity  located  at  the 
identical  position.  (The  magnitude  of  the  reflection  coefficient  varies  linearly  by  22  percent 
over  a  40  percent  frequency  bandwidth.)  Note  the  small  sideband  structure  now  evident  at 
the  base  of  the  main  impulse. 

Recall  that  we  previously  assumed  that  we  could  truncate  the  infinite  series  (1)  to  a 
finite  N-term  series  (2)  without  a  loss  in  accuracy.  In  the  spatial  domain,  this  is  the  equivalent  of 
imposing  a  finite  (spatial)  bandwidth  limitation.  That  is  to  say,  we  require  that  no  total  electrical 
path  length,  regardless  of  reflection  order,  can  be  greater  than  the  distance  lmax  (5.)  For  the 
typical  application  of  a  WR-90  waveguide-based  device,  sampled  at  256  evenly-spaced 
samples  in  p,  over  the  frequency  range  of  8  to  12  GHz,  we  compute  A I  and  lmax'  based  on  (4) 
and  (5),  to  be  2.16  and  551  inches  respectively.  In  this  particular  case  then,  the  validity  of  the 
truncation  hinges  upon  the  very  reasonable  assumption  that  all  higher-order  reflections  within 
the  device  have  dissipated  to  zero  for  path  lengths  greater  than  46  feet! 

We  mentioned  earlier  that  we  can  dramatically  reduce  data  storage  requirements  by 
exploiting  the  localization  of  spectral  energy  in  the  spatial  domain,  thereby  requiring  the 
modeling  and  storage  of  only  a  "few”  dominant  Fourier  coefficients  with  a  minimal  loss  in 
overall  modeling  accuracy.  This  economy  is  illustrated  in  Figure  13  where  the  concentration  of 
the  majority  of  spectral  energy  (for  the  two-stub  tuner)  in  a  minority  of  Fourier  coefficients  is 
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Ak  =  8  (k-i) 


FIGURE  9.  INTEGER-POSITIONED  DISCONTINUITY, 
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FIGURE  11.  FREQUENCY-INDEPENDENT  DISCONTINUITY, 
SPATIAL  DOMAIN 


FIGURE  12.  UNEARLY  FREQUENCY-DEPENDENT 
DISCONTINUITY,  SPATIAL  DOMAIN 


clearly  evident.  Figure  14  illustrates  the  loss  >n  modeling  accuracy  (in  the  frequency  domain) 
for  the  real  part  of  the  input  reflection  coefficient  for  the  two-stub  tuner  with  100  percent  (256 
terms),  95  percent  (20  terms),  and  90  percent  (13  terms)  of  the  total  spectral  energy  retained. 
As  one  might  expect,  degradation  in  modeling  accuracy  is  most  evident  at  the  inflection 
points,  where  low-pass  filtering  effects  begin  to  emerge. 

To  review  then,  referring  to  Figure  4,  data  reduction  in  the  spatial  domain  involves  the 
application  of  both  analog  and  digital  signal  processing  techniques  to  the  measured  scattering 
coefficient  data  so  as  to  reduce  the  data  to  a  compact  frequency-independent  form  (Fourier 
coefficients.)  These  coefficients  (one  set  for  each  scattering  parameter)  are  then  passed  to 
the  next  processing  stage  where  the  Fourier  coefficients  are  parametrically  modeled  in  terms 
of  their  physical  variable  dependencies. 

2.3  Physical  Variable  Dependence  Modstirg 

Figure  15  illustrates  the  overall  procedure  for  the  parametric  modeling  of  the  phyUuai 
variable  dependence  for  each  of  the  retained  Fourier  coefficients. 

For  the  case  of  the  two-stub  tuner  (pictured  in  Figure  17)  each  Fourier  coefficient  is 
dependent  upon  two  physical  variables  XI  and  X2  (the  "near"  and  "far"  stub  depths.)  Figure 
16  illustrates  the  magnitude  and  phase  "surfaces"  associated  with  the  k  =  6  Fourier  coefficient, 
corresponding  to  the  spatial  location  of  the  "far"  X2  stub.  Note  the  strong  exponential 
dependence  of  the  coefficient  magnitude  on  the  X2  variable  and  the  overall  smoothness  of 
the  surface.  In  contrast,  note  the  highly  convoluted  surface  associated  with  the  phase 
behavior  of  the  k  =  6  Fourier  coefficient.  At  first  inspection,  this  behavior  seems  peculiar  in 
light  of  the  smooth  magnitude  characteristic.  However,  noting  that  phase  (9)  is  normally 
expressed  MODULO  (9,2k),  we  recognize  that  some  of  this  ill-behavior  may  be  due  to 
"retraces"  associated  with  the  MODULO  function,  A  typical  "retrace"  in  one  dimension  is 
depicted  as  a  dashed-line  connecting  the  third  and  fourth  data  points  in  Figure  18a.  Figure 
18b  illustrates  the  same  data  after  application  of  a  DEMODULO  process  by  which  a  2rt 
constant  is  added  to  all  values  subsequent  to  the  fourth  data  point.  Clearly,  this  modified 
characteristic  is  better  suited  to  modeling.  In  fact,  application  of  simple  interpolation  between 
the  third  and  fourth  data  points  in  the  MODULO  form  would  result  in  highly  erroneous  results. 
A  two-dimensional  application  of  the  DEMODULO  process  to  the  ill-behaved  phase  data  of 
Figure  1G  is  illustrated  in  Figure  19b.  Note  the  smoothness  of  the  post-processed  surface. 


FIGURE  15.  FLOW  DIAGRAM,  MODELING  OF  PHYSICAL 
VARIABLE  DEPENDENCE  OF  FOURIER 
COEFFICIENTS 
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FIGURE  16.  PHYSICAL  VARIABLE  DEPENDENCE  OF 
FOURIER  COEFFICIENT  (A6)  FOR 
TWIN-STUB  DEVICE 


FIGURE  17. 


TWIN-STUB  DEVICE 


Given  smooth  P-dimensionai  magnitude  and  phase  (hyper-)  surfaces  for  each  Fourier 
coefficient,  it  is  a  straightforward  process  to  apply  a  P-dimensional  ieast-squares  polynomial  fit 
in  order  to  parametrically  model  the  physical  dependencies.  Weighting  factors  are  applied  to 
the  least-squares  fit  of  the  phase  (hyper-)  surface  in  order  to  achieve  local  qualities-of-fit  that 
are  proportional  to  the  Fourier  coefficient  magnitudes.  The  orders  of  these  polynomials  are 
variable  (within  the  constraints  of  the  number  of  measured  physical  variable  permutations),  and 
are  adjusted  in  order  to  maximize  modeling  accuracy  while  minimizing  the  numbei  of 
polynomial  coefficients  which  mus  be  ultimately  stored. 

2.4  Formatting  and  Storage  of  the  Scattering  Library 

With  the  parametric  modeling  process  now  complete,  the  polynomial  coefficients  for 
the  magnitude  and  phase  for  each  Fourier  coefficient  associated  with  each  scattering 
parameter  for  a  given  device  are  formatted  and  stored  in  the  scattering  library.  This 
condensed  representation  is  accompanied  by  the  device  descriptors  and  parameters 
necessary  in  order  to  reconstitute  the  scattering  parameters  at  specific  frequencies  for  specific 
physical  variable  settings  Note  that,  besides  "real"  measured  components,  fictitious  devices 
based  on  theoretical  development  can  be  synthesized,  modeled,  and  stored  as  entries  in  the 
scattering  parameter  'library. 

3.  The  Circuit  Analysis  "Engine" 

In  a  manner  similar  to  ether  software  packages,  a  schematic  representation  of  an 
arbitrary  interconnection  of  multi-port  devices  is  based  on  a  nooal  format.  Each  internal  node 
is  common  to  exactly  two  devices  with  each  external  node  common  to  only  one. 
Interconnecting  transmission  line  sections  are  specified  as  separate  two-port  devices.  (Figure 
25  illustrates  this  process.} 

Expressions  for  each  of  the  two  opposite  traveling  waves  at  each  node  can  be  written 
in  terms  of  the  scattering  coefficients  and  the  incoming  wavc-s  for  each  of  the  two  adjacent 
devices.  For  M  internal  nodes,  this  results  In  2M  equations  in  2M  unknowns  (internal  wave 


coefficients.)  We  can  therefore  solve  for  ail  internal  wave  coefficients  leaving  only  the 
scattering  parameters  for  the  N-pcrt  composite  device  defined  by  the  N  external  nodes.  This 
process  is  performed  in  practice  through  the  construction  of  a  square  (2M+N)  x  (2M+N)  matrix, 


which  combines  expressions  for  the  outgoing  waves  at  each  of  the  external  nodes  in  terms  of 
the  internal  wave  coefficients  and  the  external  incoming  waves  together  with  the 
homogeneous  expressions  for  the  internal  wave  coefficients  in  terms  of  the  external  incoming 
waves.  This  allows  for  the  direct  computation  of  the  external  scattering  parameters  by 
partitioning  the  matrix  into  four  parts  and  solving  for  the  external  outgoing  waves  in  terms  of 
the  external  incoming  waves  through  the  application  of  simple  matrix  algebra  and  a  single 
2Mx2M  inversion. 

Figure  20  illustrates  the  process  performed  by  the  circuit  analysis  engine.  Special 
care  has  been  taken  in  optimizing  the  numerical  efficiency  of  this  routine  because  of  its 
repetitive  nature  and  its  dominance  in  the  overall  execution  time  for  the  unified  analysis 
package.  Part  of  this  optimization  will  ultimately  include  the  exploitation  of  the  sparse  nature  of 
the  partitioned  matrix,  thereby  allowing  for  the  decimation  of  the  large  matrix  equation  into 
many  smaller  parts. 

4.  The  User-Specified  Device  Format 

An  obvious  prerequisite  to  the  analysis  of  a  given  circuit  is  the  specification  of  the 
properties  and  interconnections  of  the  constituent  devices.  This  is  accomplished  via  a  user- 
specified  input  file  in  which  the  node  connections  and  physical  variable  settings  for  each 
device  are  declared,  one  line  per  device.  A  typical  input  file  line  specifying  a  3-port  H-plane 
reactive  tee  is  illustrated  in  Figure  21 . 

A  device  designation  "DVDR0107"  appears  at  the  beginning  of  the  line  and  serves  as 
a  "call  number"  for  locating  the  previously  stored  and  condensed  scattering  parameters  for  this 
device  within  the  scattering  library.  The  three  integers  "15",  "11",  and  "32"  specify  the  circuit 
nodes  to  which  the  first,  second,  and  third  ports  of  the  device  are  attached.  The  following 
three  real  numbers  "120",  "3",  and  "2"  describe  the  mean  value,  absolute  tolerance,  and 
standard  deviation  for  the  physical  variable  corresponding  to  the  offset  of  the  septum  internal 
to  the  tee.  The  other  two  triplets  dictate  the  same  information  for  the  variables  associated  with 
the  septum  and  iris  depths  for  the  device.  Alternatively,  "global”  variables  can  be  specified  in 
place  of  physical  constants  in  order  to  expedite  the  adjustment  of  settings  and  tolerances  for 
similar  devices. 
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FIGURE  20.  FLOW  DIAGRAM,  CIRCUIT  ANALYSIS  ENGINE 
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FIGURE  22.  FLOW  OIAGRAM,  EXECUTIVE  PROGRAM 
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The  statistical  information  for  each  variable  describes  a  probability  distribution 
function,  truncated  at  the  tolerance  extremes.  These  specifications  are  utilized  when 
performing  statistical  studies  of  various  composite  circuit  variables  in  terms  of  the  tolerances  of 
its  constituent  devices. 

5.0  The  Executive  Program 

The  Executive  Program  is  responsible  for  the  overall  management  of  the  various 
unified  analysis  operation.,  jj'id  is  illustrated  iri  fiuw  aiagram  form  in  Figure  22.  (Refer  to  Figure 
1  for  the  location  of  the  Executive  Program  within  the  overall  unified  analysis  approach.)  Note 
that  most  of  the  features  discussed  in  this  section  are  currently  implemented  in  the  executive 
program,  while  some  remain  in  various  stages  of  ongoing  development. 

Operation  commences  with  the  reading  of  the  user-specified  schematic  and  device 
descriptions  embodied  in  the  input  file  together  with  user-specified  job  control  information 
dictating  the  desired  functions  and  options  to  be  performed.  The  input  data  set  is  checked  for 
consistency  in  the  schematic  and  device  specifications.  Verifications  include:  Is  each  node 
specified  exactly  once  (external)  or  twice  (internal)?  Are  the  number  of  nodes  in  the  node 
connection  list  equal  to  the  number  of  ports  for  the  device?  Is  the  device  model  valid  over  the 
entire  physical  variable  and  frequency  ranges  of  interest?  Are  all  port  interfaces  continuous? 

Following  error  checking,  the  user-specified  schematic  is  optimized  by  identifying 
equivalent  sections  of  the  circuit  and  therefore  eliminating  redundant  analysis  operations. 
Additionally,  the  schematic  is  decimated  into  subsections  which  are  analyzed  separately  in 
order  to  isolate  statistically  variable  devices  and  to  minimize  matrix  sizes,  thereby  enhancing 
numerical  efficiency.  Next,  device  variables  are  perturbed  based  on  their  individual  probability 
distribution  functions  (statistical  computations)  or  in  specified  differential  increments  (gradient 
computations.) 

With  all  device  variables  now  defined,  the  individual  Fourier  coefficients  associated 
with  each  scattering  parameter  for  each  device  are  reconstructed  through  the  expansion  of 
the  polynomials  whose  coefficients  are  stored  in  the  scattering  library.  Note  that  these 
reconstructed  Fourier  coefficients  have  no  frequency  dependence  and  are  therefore  valid  for 
the  entire  frequency  bandwidth  of  interest. 


The  actual  scattering  parameters  for  each  device  are  now  reconstructed  in  the 
frequency  domain  through  the  application  of  an  inverse  FFT  to  the  Fourier  coefficients  and 
the  utilization  of  an  interpolation  formula  as  dictated  by  the  Sampling  Theorem.  Interpolation 
begins  with  the  lowest  frequency  and  continues  incrementally  across  the  bandwidth  of 
interest  as  specified  by  the  user. 

As  another  optimizing  feature,  simple  transmission  line  lengths  are  "absorbed"  into 
adjacent  devices  by  altering  the  appropriate  scattering  parameter  phase  values  in  order  to 
simulate  the  presence  of  the  lines.  This  has  the  effect  of  reducing  the  number  of  required 
internal  nodes  by  a  factor  of  two  and  therefore  improves  computational  efficiency. 

Once  fully  optimized,  the  composite  scattering  parameters  for  the  (sub-)  circuit  are 
computed  via  the  circuit  analysis  engine  and  then  temporarily  stored.  This  process  of 
scattering  parameter  reconstruction,  transmission  line  absorption,  and  (sub-)  circuit  analysis  is 
repeated  for  each  frequency  of  interest. 

It  is  at  this  point  that  any  required  auxiliary  operations  are  performed  on  the  computed 
external  scattering  parameter  data.  These  operations  include  the  calculation  of  related 
network  parameters  (VSWR,  insertion  loss  and  phase,  transmission  loss  and  phase,  isolation) 
based  on  Smm  and  Smn  coefficients;  the  computation  of  far-field  antenna  characteristics 
(patterns,  gain,  RMS  sidelobe  levels)  based  on  aperture  Smn  coefficients,  lattice  spacings, 
and  element  factors;  and  the  prediction  of  "pre-aperture"  RCS  contributions  based  on 
aperture  Smm  coefficients  and  lattice  spacings.  These  characteristics  are  stored  and 
accumulated  for  all  required  statistical/gradient  perturbations  of  the  physical  variables. 

Decimated  (sub-)  circuits  are  now  reassembled  and,  based  on  the  accumulated 
results  for  the  various  perturbation  cases,  statistical  parameters  (mean  and  standard 
deviations)  and  gradients  are  computed  for  network  and  antenna  characteristics  and  then 
output  in  tabular  and/or  graphical  form  as  requested  by  the  user.  Plots  of  mean  values  and 
standard  deviation  "windows"  for  individual  characteristics  as  functions  of  frequency  and/or 
angle  are  available.  This  allows  for  the  informed  adjustment  of  the  mean  values  and 
tolerances  of  specific  physical  variables  in  order  to  meet  specified  network  and  antenna 
requirements  with  predictable  levels  of  certainty.  Gradient  information  is  used  for  optimization 
of  specific  physical  variables  subject  to  user-specified  weighting  and  penalty  functions. 
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6. 


A  Typical  Application:  Waveguide  Corporate  Feed  Network 


Figure  23  illustrates  an  eight-port  Corporate  Feed  Subassembly  comprised  of  six 
reactive  H-plane  power  dividers,  five  H-plane  bends,  and  seven  E-plane  bends,  plus 
interconnecting  waveguide  sections. 

A  generic  H-plane  power  divider  was  fabricated  as  pictured  in  Figure  24  and 
characterized  in  terms  of  its  scattering  parameters  over  a  range  of  frequencies  and  physical 
variable  combinations  (septum  offset,  septum  length,  and  iris  depth.)  Likewise,  the  two-port 
scattering  parameters  for  the  H-plane  and  E-plane  waveguide  bends  were  measured  across  a 
broad  frequency  bandwidth.  Given  the  known  physical  variable  settings  for  each  of  the  six 
reactive  H-plane  power  dividers  comprising  the  feed  network,  the  scattering  parameters  for 
each  were  reconstructed  for  all  frequencies  of  interest.  Figure  25  illustrates  the  node  and 
device  definitions  for  the  network  along  with  the  (nonstatistical)  schematic  representation  and 
the  specification  of  interconnecting  waveguide  section  lengths. 

The  unified  analysis  approach  (in  its  developmental  form)  was  applied  to  the  network 
and  Figures  26a  and  26b  illustrate  the  measured  and  predicted  insertion  losses  between  the 
input  and  each  of  the  seven  output  ports  across  the  entire  bandwidth  of  interest.  Note  the 
faithful  prediction  of  both  the  macroscopic  and  microscopic  structure  of  the  measured  data. 
Similarly,  plots  of  insertion  phase  and  input  VSWR  show  excellent  correlation  between 
predicted  and  measured  results. 

7.  Progress-to-Date  and  Future  Plans 

As  mentioned  in  the  abstract  and  introduction,  the  development  of  the  unified 
analysis  approach  is  an  ongoing  effort.  To  date,  we  have  proven  the  viability  of  our  modeling 
and  analysis  techniques  both  in  theory  and  in  direct  application  via  "proof-of-design"  versions 
of  the  requisite  software.  Direct  comparison  between  measured  and  predicted  results  for 
typical  microwave  subassemblies  have  yielded  excellent  correlations. 

We  continue  to  augment  the  scattering  parameter  library  with  the  characterizations  of 
common  devices  so  as  to  expand  the  applicability  of  our  analysis  approach  to  a  growing  range 
of  subassemblies  and  assemblies.  The  statistical  capabilities  of  the  analysis  routine  are 
presently  under  test  as  are  several  of  the  aforementioned  optimization  schemes  for  boosting 
numerical  efficiency. 
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CIRCUIT  DIAGRAM 


NODE  MATRIX 


LENGTH  MATRIX 


4.860  1.623 
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FIGURE  25.  SCHEMATIC  AND  NODAL  REPRESENTATION,  UPPER  FEED  ASSEMBLY 


FIGURE  26a. 


UPPER  FEED  ASSEMBLY  INSERTION-LOSS, 
MEASURED 


FIGURE  26b.  UPPER  FEED  ASSEMBLY  INSERTION-LOSS 
PREDICTED 


Auxiliary  capabilities  such  as  far-field  patterns  and  RCS  characterizations  are  not  as  yet 
available.  However,  implementation  of  these  features  is  expected  to  be  a  straightforward 
process  given  the  self  (Smm)  and  coupling  (S-|n)  scattering  coefficients  for  the  entire  antenna 
assembly  with  the  input  ports  and  radiating  elements  expressed  as  external  nodes. 

As  the  scattering  parameter  library  continues  to  grow  and  our  analysis  package 
becomes  further  refined,  we  will  apply  our  unified  approach  to  larger  and  larger  portions  of 
entire  antenna  systems,  with  tho  ultima’s  goal  of  modeling  their  complete  behavior,  replete 
from  sum-port  to  aperture-plate,  while  realizing  new  levels  of  analysis  accuracy  and  efficiency. 
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ABSTRACT 


The  usual  approach  in  solving  Maxwell’s  equations  with  given 
current  distribution  is  by  means  of  finding  the  vector  potential 
function  due  to  the  current  source.  This  approach  is  generally  very 
involved,  for  it  contains  vector  differential  operations  in  addition 
to  the  integration.  In  this  article,  useful  formulas  for  directly 
finding  electric  and  magnetic  fields  are  derived,  where  the 
complicated  vector  differential  operations  are  replaced  by  the  easy 
vector  algebraic  operations.  The  resulting  solutions  are  certainly 
exact  with  no  approx imat ion,  and  can  be  explicitly  written  in  matrix 
formulation  such  that  both  the  current  source  and  the  resulting 
fields  may  be  conveniently  expressed  in  terms  of  various  coordinate 
systems.  If  the  radiation  fields  in  the  near  or  far  region  are  of 
primary  interest,  the  formulation  becomes  simply  the  matrix 
mu 1 t i pi icat ion,  whenever  the  radiation  vector  is  determined. 
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SOLVING  MAXWELL  EQUATIONS  BY  MATRIX  FORMULATION 


INTRODUCTION 

In  solving  Maxwell’s  equations  with  a  given  current  distribut ion, 
the  usual  approach  Cl-51  is  to  find  the  vector  potential  due  to  the 
current  source.  The  desired  solutions  can  then  be  determined  by 
t5vng  vector  differential  operations  upon  the  vector  potential.  A 
vector  differential  operation  is  an  operation  which  contains  the 
vector  differential  operator  V,  such  as  gradient,  divergence,  curl, 
and  La p 1 ac i an.  By  introducing  the  free-space  Dyadic  Green’s  function 

Cl  I,  one  can  bypass  the  vector  potential  function,  and  thereby 
establish  a  direct  connection  between  the  current  source  and  the 
resulting  electromagnet ic  fields.  However,  the  solutions  of 
Maxwell’s  equations  are  still  quite  involved,  because  of  the  vector 
differential  operations  in  addition  to  the  integration.  In  numerical 
computation,  the  evaluation  is  usually  much  easier  for  the 
integration  than  the  di f f erent i at i on,  since  integration  can  always  be 
approximated  using  summations. 
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In  this  article  useful  formulas  for  finding  electromagnet  ic 
fields  are  csrivec.  These  formulas  replace  the  complicated  vector 
differential  operations  with  easy  vect or  algebraic  operations.  0 
vector  algebraic  operation  is  an  operation  which  contains  only  the 
dot  product  and  cross  product  of  vectors.  The  resulting  field 
solutions  are  exact  with  no  approx imat ion.  Furthermore,  the  field 
solutions  can  be  explicitly  written  in  matrix  formulation  so  that 
both  the  current  source  and  the  resulting  fields  can  be  expressed  in 
terms  of  various  coordinate  systems. 

If  the  radiating  fields  C4,5D,  either  in  the  near  field  (Fresnel 
region)  or  in  the  far  field  (Fraunhofer  region),  are  of  primary 
interest,  the  formulation  can  be  greatly  simplified.  The 
electromagnet ic  fields  are  easily  derived  by  simple  matrix 
rnul  t  i  p  1  icat  i on  whenever  the  radiatiori  vector  L4]  is  calculated. 

Several  examples  appearing  in  electromagnet ic  theory  will  be 
given  to  demonstrate  the  merit  of  the  technique  presented. 
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FORMULATION 


It  is  well-known  Cl-51  that  the  solutions  of  Maxwell  eouations 
in  time  harmonics  form 
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where  A  is  called  a  vector  potential,  related  to  the  given  current 
distribution  J  over  a  volume  V, 
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G y<r ) g ( I r  -  r  •  I)  is  the  well-known  free— space  Dyadic  Green’ s  operator. 
The  electrical  field  and  magnetic  field  may  therefore  be  written  as 
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Evaluation  of  (2)  or  (6)  is  in  general  not  very  easy,  since  it 
involves  both  vector  differential  operations  and  integrat ions.  It  is 
shown  in  the  Appendix  that  the  complicated  vector  differential 
operations  in  (6)  can  be  eliminated  and  replaced  by  easy  vector 
algebraic  operations-  The  results  are 
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where  G(R)  and  K ( R )  are  vector  algebraic  operators,  containing 
only  the  dot  product  and  cross  product  of  vectors. 
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K  ( R) 
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1 

kR 


R  x 
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and  R  is  a  radial  vector,  independent  of  the  coordinate  system, 
directed  along  the  line  of  sight  from  the  source  position  r'  to  the 
observation  position  r,  as  illustrated  in  Figure  1(a). 

R  =  r  _  r ' ,  R  =  I  r  -  r '  I  ,  R  =  — L  (9) 

R 


Both  (6)  and  (7)  are  exact  solutions  of  Maxwell’s  equations  (1>, 
directly  related  to  the  current  source.  The  evaluation  by  (7)  is  more 
manageable  now  compared  to  (6),  since  it  involves  only  integration. 
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To  demonstrate  this,  we  may  convert  the  vector  expression  into  matrix 
expression  by  using  the  technique  of  matrix  formulation  of  vector 
operations  C6D.  Equations  (7)  may  then  be  written  as, 
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(10b) 


V 

where  the  spherical  coordinates  s  (0R,  R)  is  defined  such 

that  R  is  a  radial  vector  from  r'  to  r.  Pnd  a  vector,  such  as 
(E(rP))P,  is  a  function  of  variables  in  the  q -coord inaites  and  has 
three  components  in  the  p-coord i nates.  ft Iso  a  coordinate 
transformat  ion  matrix  relating  the  u— coord inates  to  the  v-coordinates 
is  denoted  by  CVTUD,  which  can  readily  be  fou.id  by  the  technique 
of  a  coordinate  transformation  flowchart  C7D.  fts  a  practical 
application,  the  solutions  E(r)  and  H(r)  are  usually  expressed  in 
the  spherical  coordinates  s  :  (6,  0,  r) ,  and  the  source  J(r' )  may 

be  given  in  the  Cartesian  coordinates  c' :  (x',y',z'>.  Then  (10) 

becomes. 


_  s  f  s  S  .  ^  S  S  C'  _  .  C' 

(E  (rS)  )  -  —  j  kr^  g(R)  C  T  DC  G(R)3  C  T  D  (J(rC)>  dv'C  (11a) 


(H(rs)  ) 


-jk  j  g  ( R 


S  SR  =  -  *R  SR  c'  -  c'  c'  c' 

)  C  T  "DC  K  (R)  D  "c  T  D  (J(r'  >>  dv'C  (lib) 


Equations  (11)  are  then  written  explicitly  in  the  matrix  formulation. 
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The  solutions  of  Maxwell’s  equations  given  by  (12)  are  exact.  The 
expression  in  matrix  form  is  the  most  convenient  for  numerical 
computat ion. 

In  general,  the  radiating  field  C5D,  either  in  the  near  field  or 
in  the  far  field,  is  of  primary  interest.  If  r'  «  r  ~  X  ,  then 

R  -*■  r,  R  ■+•  r,  0r  0,  and  R  -*-  r  for  the 

amplitude  factor  and  R  r  —  r  •  r'  for  the  phase  factor.  The 

vector  potential  given  by  (3)  then  becomes 


E  (r)  =  -jkr|  G^(r)  9  N  (r> 

H  (r)  =  -jk  K^(r)  g  <r)  N  (r> 
and  from  (7), 


E  (r)  =  -jkr^  g  (r)  G  (r)  N  (r) 
H  (r)  =  -jk  g  (r)  K  (r)  N  (r> 


1  <<  X  (13a) 

r '  <<  r  — •  x  (13b) 

( 14a) 

(14b) 
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Cotl.  (14)  and  (15>  should  lead  to  the  same  results,  when  the  vector 
potential  A(r)  or  the  radiation  vector  N(r)  (for  r'  <<  r  )  is 
computed.  It  is  interesting  to  note  that  the  solution  by  (14)  is 
often  extremely  complicated,  because  it  involves  vector  differential 
operations  on  a  vector  function  A(r)  =  g(r)N(r>.  On  the  other  hand, 
the  solution  by  (15)  is  surprisingly  simple  and  straightforward, 
because  it  involves  only  matrix  multiplication.  This  can  be  seen 
clearly  from  the  following  matrix  formulation: 


From  (16)  or  (17)  we  see  that  for  r'  «  r,  the  desired  solutions 
E(r)  and  H(r)  of  Maxwell's  equations  can  be  easily  determined 
whenever  the  radiation  vector  N(r)  has  been  computed. 

To  compute  the  radiation  vector  N(r)  in  general  cases,  we  assume 
that  the  coordinate  systems  of  the  source  ~r'  and  the  field  r  are 
not  coincident,  but  are  related  through  both  origin  translation  along 
the  offset  vector  rb  *  *b'  Yb'  zb>  ant^  three  rotations  via 
Euler  angles  ($w,  €>w, 'S'w*  as  depicted  in  Figure  1(b)  The  desired 
radiation  vector  may  be  obtained  by  either  of  the  following  forms: 
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It  is  more  useful  when  we  write  (18)  explicitly  in  matrix  form. 


/'NQ(©,  ©,  r) 

© 

/ 

jKp 

N  (©,  ©,  r ) 

U 

=  e 

© 

^  N  (G,©,r)  j 

V 

' cosGcos©  cos@sin0  -sinG  ' 

-sin©  cos©  0  J 

l 

i 

sinGcos©  sinGsin©  cosG  ) 


/cos©  cos©  costji  -sin©  sinxk  — cosG  cos©  sirrfc  -sin©  cosijr  sin©  cos  © 

WWW  w  w  w  w  w  w  w  w  w 

cosG  sin©  cosifr  -*cos(f>  sin*  -cosG  sin©  sin*  +cos©  cos*  sin©  sin©  I 

w  W  W  W  W  W  W  W  W  W  W  W  ] 


-s i nG  cos* 

w  w 


sin©  si n* 

w  w 


cos© 


either 


Jx,<x,  y;  z)  ''j 

I 

Jy.(X’’  y* z) 

^  J_.<x»yi  z‘>  • 


e^ kp  dx ' dy * dz ' 


( 19a) 


cos©'  — s  i  n©'  0 

s  i  n©'  cos©'  0 

0  0  1 


J  (pi  ©,  z>  \ 

P  I 


J  ,(p’,©,  z)  e^kp  pdp'd©dz‘  (19b) 


Jz><p'i  <*>\  *>  J 


Jx.(p',©',  z) 

Jy  .(P’  ♦'»  2  ) 

JZ.<P»  Pi  z  > 


pdpd©dz 


(19c) 


,'cos©  -sin©'  0 
sin©'  cos©‘  0 
,  0  0  1 


J  ,(x  -,  yi  z ) 

P 

ja(x»  y» z> 
© 

J_.<x'i  y»  z> 


e^kp  dxdy’dz'  (19d> 


400 


where 
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In  the  case  that  there  is  no  origin  translation  and  no  axis 
rotations,  such  as  depicted  in  Figure  l^a.  then  (19)  become  simply. 
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It  is  ofteri  quite  convenient  to  include  both  an  electric  source  J 
and  magnetic  source  M  in  Maxwell’s  equations  Cl— 53: 

V  x  H  =  jw«E  +  J  (20a) 

-  9  x  E  =  j^yuH  +  M  (20b) 


The  solutions  are  then  the  superposit ion  of  the  fields  due  to  both 
current  sources. 


where 
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are  electric  and  magnetic  vector  potentials. 

If  the  radiation  field  is  of  primary  interest  (  r' <<  r  ~  X  >, 
then  the  solutions  can  be  simplified  to 
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are  radiation  vectors  due  to  the  respective  current  sources. 
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As  before  (23)  contains  only  vector  algebraic  operations  which  are 
convenient  for  numerical  computations.  We  therefore  have  the  formal 
solution  for  any  problem  consisting  of  electric  and  magnetic  currents 
in  an  unbounded  homogeneous  region.  These  formulas  include,  by 
implication,  sheets  and  filaments  of  the  electric  and  magnetic 
currents. 


EXAMPLES 

Several  examples  will  be  presented  to  demonstrate  the  merit  of 
the  technique  derived. 

1.  Hertzian  dipole  antenna 

The  current  distribution  of  an  electric  dipole  may  be  expressed 
as 
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The  corresponding  radiating  fields  are  therefore  obtained  from  (16) 
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£.  Dipole  above  the  ground  plane 

A  horizontal  dipole  is  situated  above  an  infinite  conducting 
plane.  Employing  the  image  method,  the  location  (xb,  y^,,  z^)  and 
orientation  <4>w,  ©w,  *Jfw>  are,  respect  i  vely,  (0,  0,  h)  and 
<0\  90°,  0°  )  for  the  real  element  and  (0,  0,  -h)  and  (0^  -90°,  0°  ) 

for  the  image  element.  The  current  distribution  of  both  elements  is 
assumed  to  be  that  of  the  electric  dipole  in  the  previous  example. 

By  superposit ion  and  from  (19a),  we  have 
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The  radiating  fields  are  then  followed  from  (16)  or  (17). 
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3.  Circular  loop  current 

The  current  distribution  of  a  circular  current  loop  of  radius  a 
may  be  expressed  as 
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Specializing  for  the  far  field,  where  r  >>  x » 


the  radiating  fields 


are  then  from  (17), 
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4 


Helix  antenna 


The  current  distribution  of  a  helix  antenna  as  shown  in  Fig.  £.11 
of  C33  may  be  expressed  as 
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where 

h(0)  *  (  ^  - j  U  )  J  a2  +  b2  -  kbcos© 
These  are  the  same  results  as  Eq.  2.50  -51  of  C3D. 
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5.  Circular  aperture  radiator 

The  exact  solutions  of  Maxwell’s  equations  with  given  magnetic 
current  distribution  are  similar  to  those  of  (10).  Let  the  fields  be 
in  the  Cartesian  coordinates,  and  the  source  be  a  function  of 
variables  in  the  cylindrical  coordinates  with  three  components  in  the 
Cartesian  coordinates.  Then  the  desired  solutions  in  the  matrix 
formulation  are 
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R  =  C  (x  -  p'  eosp  V  +  (y  -  p'  sinp'  )*  ■*-  (z  -  z'  >2  3  /2 

6  =  cos  *  C  (z  -  z'  )  /  R  3 

R 

p  =  tari  *C  (y  -  p'  sinp'  )  /  (x  -  p' cosp'  )  3 

R  • 


For  a  circular  aperture  of  radius  a  with  x—  polarized  magnetic 
current  distribution. 
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For  the  case  of  a  uniform  current  distribution,  we  can  evaluate  the 
integrals  in  a  closed  form  for  observation  points  along  the  axis.  To 
show  this,  let  K(p',  <f>‘  )  =  1,  and  <x,y,  z)  =  (0, 0,  z).  The  integrals 
in  z'  and  <f>‘  can  readily  evaluated.  The  results  are. 
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where  change  of  variable  by  a'-  jk(p'2  +  z':)^,  and  S  -  <z2  +  a2  )  ‘ 
Performing  the  integration  by  parts  yields. 
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These  are  the  exact  solutions  in  the  closed  for*fct  expressions  the  same 
as  those  given  by  eqs. (8)  and  (11)  of  [83. 
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CONCLUSION 


The  useful  formulas  for  solving  Maxwell's  equations  with  given 
current  distribution  are  derived.  The  results  are  exact  with  no 
approximation,  and  can  be  explicitly  written  in  matrix  formulation, 
such  that  both  the  current  source  and  result ing  fields  may  be  easily 
expressed  in  terms  of  various  coordinate  systems.  If  the  radiation 
fields  either  in  the  near  or  far  region  are  of  primary  interest,  the 
formulation  is  greatly  simplified.  The  electric  and  magnetic  fields 
can  be  simply  found  by  the  direct  matrix  multiplication  whenever  the 
radiation  vector  due  to  the  current  source  is  determined.  Some 
well-known  antenna  conf i gurat ions  are  used  to  demonstrate  the  merit 
of  the  technique  presented. 


APPENDIX 


For  reference,  some  useful  formulas  are  listed  below.  These 
formulas  are  independent  of  coodinate  system.  Also  J  is  treated  as 
a  constant  under  the  vector'  differential  operator  9. 
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where,  for  instance,  both  r 


arid  r  '  are  in  the  Cartesian  coordinate 
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The  derivation  of  (7)  from  (£)  may  therefor  be  made  by  using  some 
of  the  above  formulas,  such  as  (£7)  and  (£8). 


Then 
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This  completes  the  derivation. 
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<0)  VECTOR  OPERATORS 


1.  Vector 

2.  Dot -Product 

3.  Cross-Product 


4.  Gradient 

5.  Divergent 

6.  Curl 


7.  Divergent -Cur 1 

8.  Curl -Gradient 

9.  Dx verQent- Grad ient 

(Scalar  LaPlacian) 


10.  Gradient-Diver9ent 


11.  Curl -Curl 


12.  Grad ient -Di ver9ent 

(Vector  LaPlacian) 
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9  x  9  x 
9  2 


(A)  CARTESIAN  COORDINATES 
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<B)  CYLINDRICAL  COORDINATES; 
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(C)  SPHERICAL  COORDINATES;  s  =  (  6«  r  ) 
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